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ABSTRACT 


This  study  develops  cost  estimating  relationships  (CER)  to 
estimate  the  cost  of  supply  support  for  the  maintenance  of 
aircraft  simulators.  These  CER's  will  be  used  in  the 
Contractor  Operation  and  Maintenance  of  Simulators  (COMS) 
environment.  The  analysis  techniques  used  to  develop  and 
validate  these  relationships  include  linear  regression  (simple 
and  multiple),  analysis  of  univariate  distributions, 
simulation  modeling,  and  linear  optimization  modeling.  The 
regression  analysis  concludes  that  no  useful  CER  is  present  in 
the  cost  elements  represented  by  the  sample  of  data  gathered 
from  actual  simulator  operations.  There  is,  however,  a  useful 
CER  present  in  another  smaller  set  of  data  derived  from  the 
successful  bids  of  eight  contractors.  These  submissions  were 
for  maintenance  with  and  without  supply  support.  A  simulation 
model  was  constructed  to  provide  an  independent  cost  estimate 
for  use  with  the  maintenance  CER  developed  above.  Analysis  of 
univariate  distributions  was  used  to  transform  maintenance 
data  from  an  operational  simulator  suite  for  use  in  the 
verification  of  the  simulation  model.  Finally,  an 
optimization  model  formulation  is  recommended  for  further 
investigation  to  determine  the  best  mix  of  contractor-provided 
and  Navy-provided  spare  parts  to  complete  the  optimization  of 
Navy  expenditures  for  supply  support  for  aviation  simulators. 


V 


I .  INTRODUCTION 

A.  BACKGROUND 

In  1969  all  aviation  simulators  were  maintained  and 
operated  by  personnel  with  a  Tradevman  (TD)  rating.  The  Navy 
supply  system  was  responsible  for  providing  all  repair  parts 
and  replacement  assemblies  for  simulators.  This  system 
provided  an  adequate  number  of  operational  ready  (OR) 
simulator  hours  to  meet  the  needs  of  the  various  aviation 
communities  throughout  the  Navy. 

In  the  early  1980 's,  however,  the  Department  of  the  Navy 
dramatically  changed  the  way  aviation  simulators  were 
maintained  and  operated.  The  Defense  Officer  Personnel 
Management  Act  (DOPMA) ,  force  level  constraints  and  the  goal 
to  increase  the  size  of  the  Navy  to  600  ships  were  the  major 
factors  that  prompted  the  changes.  In  order  to  man  600  ships, 
many  new  "at  sea"  billets  would  be  needed,  and  the  Navy  end 
strength  was  already  at  the  maximum  level  permitted  by  law. 
Decisions  by  the  Congress  forced  these  additional  billets  to 
be  generated  internally.  The  additional  "at  sea"  manpower  was 
generated  by  replacing  much  of  the  military-manned  permanent 
shore  services  with  civilian  commercial  contracted  services. 

The  maintenance  and  operation  of  simulators  became  a  prime 
target  for  conversion  to  commercial  contract,  and  the  concept 
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of  Contractor  Operation  and  Maintenance  of  Simulators  (COMS) 
was  born. 

The  conversion  to  COMS  started  in  earnest  in  1984  and  was 
virtually  complete  by  1985  when  the  TD  rating  was 
disestablished.  The  original  contracts  were  for  maintenance 
and  operation  only,  with  the  Navy  continuing  to  be  responsible 
for  supply  support.  The  COMS  contracts  were  written  so  that 
the  contractor  was  not  charged  down  time  when  the  Navy  was  not 
able  to  supply  the  required  repair  parts. 

The  maintenance  of  a  complete  stock  of  repair  parts  is 
difficult  for  the  Naval  supply  system.  The  large  number  of 
parts  unique  to  simulators,  coupled  with  their  low  utilization 
rates,  places  an  undue  burden  on  a  system  which  is  not 
designed  to  handle  parts  on  a  small  scale.  The  DOD 
acquisition  regulations  further  hamper  efforts  to  have 
sufficient  spare  parts  to  support  the  COMS  contractor.  The 
requirement  for  competition,  or  waiver  for  sole  source 
acquisition,  and  the  lead  time  required  by  the  contracting 
process  itself,  frequently  cause  unusually  high  prices  and 
unacceptable  delays  in  the  arrival  of  vital  repair  parts.  The 
effect  of  these  delays  can  cause  an  excessive  number  of 
unusable  trainers  and  lead  to  decreased  readiness  for  Naval 
aviation  units. 

One  solution  to  this  dilemma  may  be  a  new  generation  COMS 
contract  that  makes  the  contractor  responsible  for  some  level 
of  supply  support.  To  aid  in  the  conduct  of  a  feasibility 
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study  to  support  the  implementation  of  this  new  COMS  contract, 
it  is  necessary  to  devise  methods  to  answer  the  following 
questions: 

1.  How  much  does  it  cost  to  provide  supply  support? 

2.  What  is  the  cost  of  a  COMS  contract? 

3.  Is  there  a  best  mix  of  contractor-supplied  or 
Navy-supplied  spare  parts? 

B.  SCOPE  AMD  LIMITATIONS 

1.  Scope 

This  research  focused  on  the  development  of  two  CERs 
that  can  be  used  to  estimate  the  cost  of  spare  parts  and  the 
cost  of  a  COMS  contract  required  to  maintain  a  simulator  for 
one  year. 

This  thesis  summarizes  the  analysis  and  presents  the 
results.  Additional  ideas  are  presented  that  outline  the 
procedures  needed  to  determine  the  best  mix  of  Navy 
furnished/contractor  furnished  spare  parts. 

2.  Limitations 

The  accuracy  of  the  results  of  this  analytical  process 
is  no  better  than  that  of  the  data  used.  Further  limitations 
exist  because  of  the  small  number  of  available  simulators  with 
useful  data  for  this  project.  Other  specific  limitations  will 
be  discussed  throughout  this  thesis  as  they  occur. 

The  intent  of  this  thesis  is  to  demonstrate  analytical 
procedures  that  lead  to  the  determination  of  a  CER.  The 
accuracy  of  the  resulting  CERs  is  not  expected  to  be 
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sufficient  to  allow  its  use  in  the  contracting  or  budgeting 
process . 
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II.  COST  ESTIMATION 


A  cost  estimate  is  a  judgment  or  opinion  regarding  the 
cost  of  an  object,  commodity,  or  service.  This  judgment  or 
opinion  may  be  arrived  at  formally  or  informally  by  a 
variety  of  methods,  all  of  which  are  based  on  the 
assumption  that  experience  is  a  reliable  guide  to  the 
future."  [Ref.l:p.l] 

Cost  estimation  techniques  are  diverse  and  range  from 
intuition  to  detailed  application  of  labor  and  engineering 
cost  standards.  The  five  techniques  most  prevalent  in  The 
Department  of  Defense  are  industrial  engineering  standards; 
rates,  factors,  and  catalog  prices;  cost  estimating 
relationships;  specific  analogies;  and  expert  opinion.  There 
are  many  more  techniques  in  the  civilian  sector.  The 
determination  of  the  technique  to  be  applied  is  dependent  on 
the  data  available,  what  the  estimate  is  to  be  used  for,  the 
time  available  to  do  the  study,  and  other  factors  that  affect 
the  direction  of  the  study.  For  the  purposes  of  this  thesis, 
cost  estimating  relationships  and  industrial  engineering 
standards  will  be  the  only  two  methods  explored. 

A.  COST  ESTIMATING  RELAXIONSHIFS 

Cost  estimating  relationships  (CERs)  are  equations  derived 
from  statistical  analysis  that  are  used  to  estimate  the  cost 
of  a  product  based  on  some  characteristic  of  that  product. 
This  method  is  preferred  because  of  its  ease  of  application 
and  accuracy.  Its  use  requires  the  availability  of  high 
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quality  data  and  the  existence  of  relationships  among  elements 
of  that  data  base.  A  CER  can  be  generated  using  regression 
analysis  processes.  The  details  of  that  process  are  discussed 
in  the  next  chapter. 

B.  INDUSTRIAL  ENGINEERING  STANDARDS 

On  some  occassions  sufficient  data  are  not  available,  or 

a  useful  CER  can  not  be  derived  from  the  data,  or  a  change  in 

the  process  or  materials  makes  historical  data  unusable.  In 

these  instances  industrial  engineering  standards  must  be  used. 

Estimating  by  industrial  engineering  procedures  can  be 
broadly  defined  as  an  examination  of  separate  segments  of 
work  at  a  low  level  of  detail  and  a  synthesis  of  the  many 
detailed  estimates  into  a  total.  [Ref.  l:p.2] 

In  this  estimation  scheme  the  process  of  production  is 
broken  down  into  its  smallest  elements.  The  cost  of  labor  and 
material  for  each  of  these  elements  of  production  is  then 
computed  and  totaled.  The  results  are  then  multiplied  by 
factors,  determined  by  management,  to  estimate  overhead  and 
other  indirect  costs.  The  summation  of  these  direct  and 
indirect  costs  equals  the  total  cost  estimate  for  the  product. 

C .  PROBLEMS 

As  is  the  case  with  all  estimation  processes,  there  are 
some  inherent  problems  in  each  of  these  procedures.  In  the 
case  of  cost  estimating  relationships,  there  are  two 
limitations.  The  first  is  the  inherent  uncertainty  always 
associated  with  the  use  of  statistics.  The  second  is  the 
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uncertainty  of  the  application  of  a  particular  CER  to  a 
certain  situation. 

The  industrial  engineering  standards  technique  has  three 
drawbacks.  First,  it  takes  a  large  number  of  personnel  with 
a  detailed  knowledge  of  the  production  process  to  compute  the 
production  element  costs.  Second,  small  mistakes  in  the 
production  element  costs  can  be  magnified  when  mulciplied  by 
the  predetermined  factors  in  the  computation  of  indirect 
costs.  Third,  the  whole  cost  frequently  turns  out  to  be 
greater  than  the  sum  of  all  the  production  element  costs. 

In  considering  the  benefit  to  be  derived  from  a  cost 
estimation,  one  must  always  keep  in  mind  that  the  estimations 
must  be  reasonable  and  structurally  sound.  The  usefulness  of 
the  cost  estimation  may  also  be  affected  by  the  analyst 
himself.  The  lack  of  personal  knowledge  or  familiarity  with 
the  systems  being  considered  or  a  personal  bias  interjected 
into  the  analysis  by  the  analyst  can  detract  from  the 
credibility  of  the  estimate. 
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III.  METHODS  OF  ANALYSIS 

A.  DETERMINING  COST  ESTIMATING  RELATIONSHIPS 
1.  Simple  Regression  Analysis 

Simple  regression  analysis  is  defined  and  explained  in 
books  that  discuss  data  analysis.  In  Bachelder's  An 
Introduction  to  Equipment  Cost  Estimating  [Ref.  1],  simple 
linear  regression  is  defined  as  a  simple  two-variable  model 
which  describes  the  linear  relationship  between  the  two 
variables  by  the  equation 

y=a+bx+e  (3.1) 

where  y  represents  the  dependent  variable  (cost  in  this  case) , 
X  represents  the  independent  variable  (any  of  the  other  data 
elements) ,  and  e  is  a  random  error  or  fluctuation.  The 
constants  a  and  b  represent  the  Y-axis  intercept  and  the  slope 
of  the  line,  respectively.  Equation  (3.1)  represents  the 
linear  relationship  between  x  and  y,  and  it  can  be  used  to 
estimate  y  if  an  accurate  value  is  available  for  x.  The 
parameters  a  and  b  are  estimated  from  data  pairs  (Xj, 
yj)...(x^,  y^)  using  well  established  analytical  methods.  The 
equations  for  these  estimates  A  and  B  for  a  and  b, 
respectively,  are  as  follows: 
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A=y~Bx 


Tix^-x)  (Vj-y) 
’Lix^-x)^ 


when 


n 


5^y2 .  y=  —^y- 


n 


The  STSC  program  Statgraphics  was  used  in  this  thesis  to 
perform  the  statistical  analyses.  This  program  uses  the 
ordinary  least  squares  method  to  estimate  the  values  of  a  and 
b,  and  will  perform  the  calculations  need  to  make  statistical 
inferences  about  the  parameters  a  and  b  and  the  line  a  +  bx. 


a.  Standard  Error  of  the  Regression  ?SE) 

The  standard  error  of  the  regression  is  a  value 
that  represents  the  standard  deviation  of  the  estimated  values 
of  the  error,  e.  The  equation  for  SE  in  terms  of  the  data 
pairs  (Xj,  y,)  is 

SE=^'L{y^-  (A^Bx^))^ 


where  A  and  B  are  the  least  squares  estimate  of  a  and  b, 
respectively. 
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b.  T-value 


defined  as 


The  t-value  associated  with  the  parameter  b  is 


SE/'LiXi-x)'^ 

It  is  used  to  test  the  significance  of  the  relationship 
between  x  and  y.  This  is  done  by  performing  a  statistical 
test  of  the  hypothesis  that  the  slope  coefficient,  b,  is  equal 
to  0.  That  is,  it  is  decided  that  h  ^  0  when  the  value  of 
the  t-statistic  in  equation  3.2  is  larger  in  absolute  value 
than  the  value  of  a  number,  t^  taken  from  standard 
student-t  tables.  The  value  of  a  is  the  level  of  significance 
for  this  test  of  hypothsis. 

c.  R-scmared  Value 

The  R-squared  value  is  a  measure  of  how  well  the 
regression  equation  predicts  the  value  of  y.  It  is  frequently 
called  the  coefficient  of  determination,  and  equals  one  if  all 
of  the  observed  values  are  on  the  predicted  line  of  the 
estimated  model,  A  +  Bx. 

2.  Multiple  Regression 

If  the  regression  of  one  variable  on  another  fails  to 
produce  an  adequate  CER,  then  a  technique  called  multiple 
regression  can  be  used.  This  method  is  similar  to  simple 
regression  except  it  has  more  than  one  independent  variable. 
In  the  resulting  representative  equation 
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y=a+bx+cz+e,  (3.3) 
y  is  still  the  dependent  variable,  x  and  z  are  independent 
variables,  e  is  the  random  error  or  fluctuation,  a  is  the 
Y-axis  intercept,  and  the  coefficients  b  and  c  represent  the 
change  in  y  for  each  change  of  one  unit  of  their  respective 
variables  (x  and  z) . 

It  is  important  to  remember  that  each  added  variable 
reduces  the  accuracy  of  the  statistical  inferences  made  about 
the  paramters  a,  b,  c  and  about  the  line  itself.  This  becomes 
critical  when  working  with  small  data  sets  as  it  can  decrease 
the  credibility  of  the  analysis. 

The  computed  t-value  and  R-squared  numbers  are  used 
for  the  same  purposes  as  in  simple  regression. 

a.  Collinearity 

Collinearity  exists  when  two  of  the  variables  in 
a  multiple  regression  model  have,  between  themselves,  a  linear 
relationship.  An  indication  that  this  condition  exists  is  the 
presence  of  a  high  R-squared  value  while  one  or  more  of  the 
variables  has  a  low  t-value. 

B.  SIMULATION  MODELING  TECHNIQUES 

If  statistical  analysis  fails,  another  form  of  analysis 
known  as  simulation  modeling  may  be  used.  Simulation  modeling 
is  the  second  technique  to  be  utilized  in  this  thesis. 

One  dictionary  definition  of  the  verb  simulate  is  as 
follows:  "to  assume  or  have  the  appearance  or  characteristics 
of."  [Ref.  2] 
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This  would  appear  to  be  a  reasonably  accurate  description  of 
the  technique.  The  key  to  a  successful  simulation,  however, 
is  not  just  the  "appearance  or  characteristics,"  but  rather 
the  operation  or  behavior  of  the  simulation  model. 

A  useful  working  definition  of  simulation  is  the 
following: 

The  process  of  designing  a  mathematical  or  logical  model 
of  a  real  system  and  then  conducting  computer-based 
experiments  with  the  model  to  describe,  explain,  and 
predict  the  behavior  of  the  real  system.  [Ref.  3:p.  5] 

Utilization  of  this  second  technique  requires  the 
construction  of  a  model  of  the  aircraft  simulators  that 
simulates  the  maintenance  process.  This  model  utilizes 
param.eters,  unique  to  the  performance  of  the  system  being 
investigated,  to  estimate  the  number  of  man-years  required  to 
maintain  the  system.  This  information  is  then  used  to 
determine  the  cost  of  a  COMS  contract  for  a  year. 

To  validate  this  model,  it  is  necessary  to  analyze 
historical  data  from  a  known  system  and  use  the  analysis 
results  to  run  the  simulation.  The  results  of  the  simulation 
are  compared  to  reality.  If  the  two  are  statistically 
equivalent,  the  model  is  validated  and  usable. 

1.  Analysis  of  Maintenance  Data 

The  analysis  of  the  maintenance  data  is  fairly 
straightforward.  First,  it  is  necessary  to  obtain  a  data  base 
that  contains  the  time  each  failure  occurred  and  the 
maintenance  time  it  took  to  repair  each  discrepancy.  The  time 
of  failure  data  is  used  to  determine  the  inter-arrival  times 
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(time  between  failures) .  The  inter-arrival  times  and  repair 
times  are  then  analyzed  to  determine  their  probability 
distributions  and  the  parameters  of  those  distributions. 
These  parameters  are  fed  into  the  model  which  then  simulates 
the  maintenance  process. 

a.  Hi 

The  first  step  in  analysis  is  to  display  a 
histogram  of  the  data.  A  histogram  is  a  graph  that  shows  the 
number  of  occurrences  of  data  points  within  a  bin  (range  of 
values) .  The  histogram  will  indicate  the  location  (mean)  and 
variance  of  the  probability  distribution  associated  with  the 
data.  It  may  also  give  additional  indications  about  this 
probability  distribution. 

b.  Probability  Plot 

The  probability  plot  is  the  primary  means  of 
determining  the  distribution  reflected  in  the  data  set.  This 
plot  compares  the  quantile  of  a  data  set  with  the  theoretical 
quantile  of  a  selected  distribution.  If  the  two  plots  are 
reasonably  close,  the  probability  distribution  that  generated 
the  sample  is  said  to  be  the  same,  statistically,  as  the 
theoretical  distribution. 
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c.  Kolmogorov  -  Smirnov  Test 


The  Kolmogorov  -  Smirnov  test  is  used  to  determine 
the  goodness-of-f it  of  a  sample  distribution  to  a  theoretical 
distribution.  It  does  this  by  comparing  the  absolute  maximum 
distance  between  the  sample  cumulative  distribution  function 
and  the  theoretical  cumulative  distribution  function.  The 
computed  significance  level  is  used  to  determine  if  the  sample 
distribution  is  the  same  as  the  theoretical  distribution. 
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IV.  COST  ESTIMATION  MODELS 


A.  COST  DATA 

The  data  used  in  this  thesis  were  made  available  by  Fleet 
Aviation  Specialized  Operational  Training  Group,  Atlantic 
Fleet  (FASO) .  The  data  represent  the  annual  costs  incurred  in 
the  acquisition  and  operation  of  forty-six  separate  aviation 
simulators  for  FY-87.  The  data  elements  provided  include  the 
cost  of  device  acquisition  (ACQCOST) ,  number  of  hours  the 
simulator  was  programmed  for  use  (PGMHRS) ,  COMS  contract  cost 
(COMCOST) ,  utilities  cost  (UTILCOST) ,  spare  parts  cost 
(SUPPCOST) ,  total  cost  (TOTLCOST) ,  and  the  cost  per  hour  to 
operate  the  simulator  (HRLYCOST) .  TOTLCOST  is  the  sum  of 
COMCOST,  UTILCOST,  SUPPCOST,  and  overhead  costs  for  contract 
compliance  inspection  personnel.  HRLYCOST  is  TOTLCOST  divided 
by  PGMHRS. 

Additional  data  on  the  cost  of  eight  five-year  COMS 
contracts  bid  with  and  without  supply  support  were  provided  by 
Naval  Training  Systems  Center  (NTSC) ,  Orlando,  Florida. 

The  first  step  in  the  analysis  of  this  data  was  to  see  if 
a  CER  existed  that  would  estimate  the  cost  of  supply  support. 
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B.  DEVELOPMENT  OF  CER  USING  FASO  DATA 


1.  Simple  Regression 

Using  the  cost  data  supplied  by  FASO,  a  simple 
regression  of  SUPPCOST  against  all  the  other  data  elements 
demonstrated  the  absence  of  a  useful  CER  with  any  of  them. 
For  the  purposes  of  comparison,  the  t-statistic  for  45  degrees 
of  freedom  and  level  of  significance  a=.05  is  2.414.  A 
statistical  summary  of  these  regressions  is  provided  in  TABLE 
4.1.  Note  that  the  t-value  of  the  slope  in  all  cases  but 
PGMHRS  is  larger  than  the  t-statistic  indicating  that  they  are 
all  significant,  implying  that  the  b  is  not  egual  to  0. 
However,  the  R-sguared  values  are  all  much  lower  than  that 
reguired  for  a  reasonable  explanation  of  the  dependent 
variable . 

TABLE  4.1  SIMPLE  REGRESSION  OF  INDEPENDENT  VARIABLES  AGAINST 
THE  DEPENDENT  VARIABLE  SUPPCOST 


INDEPENDENT  PARAMETER  T  R 

VARIABLE _ VALUE _ SQUARED 


UTILCOST 

SLOPE 

HRLYCOST 

SLOPE 

TOTLCOST 

SLOPE 

ACQCOST 

SLOPE 

PGMHRS 

SLOPE 

COMCOST 

SLOPE 

2.  Multiple  Regression 

With  the  failure  of 
reasonably  accurate  CER, 


4 . 88562 

35.17% 

7.68515 

57.31% 

9.61059 

67.73% 

4.94474 

35.72% 

.42723 

.41% 

3.60911 

22.84% 

simple  regression  to  produce  a 
a  backward  stepwise  multiple 
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regression  technique  was  tried.  Since  SUPPCOST  is  a  component 
of  both  TOTLCOST  and  HRLCOST,  it  would  have  a  linear 
relationship  with  both  of  those  variables.  Therefore,  in 
order  to  insure  that  collinearity  was  not  a  problem,  TOTLCOST 
and  HRLCOST  were  not  considered  during  the  multiple 
regressions . 

TABLE  4.2  shows  the  results  of  the  multiple  regression  of 
SUPPCOST  against  the  remaining  four  variables.  The  first  step 
yields  little  significance  for  any  of  the  variables,  and  the 
R-squared  value  is  very  low  at  34.21.  The  next  step 
eliminates  the  variable  with  the  lowest  t-value  (COMCOST)  and 
does  the  regression  process  again.  The  results  of  this  step 
are  not  much  better.  The  significance  of  each  of  the 
variables  is  less,  and  the  R-squared  value  does  not  change 
appreciably.  These  results  indicate  that  these  data  do  not 
exhibit  the  relationships  required  to  produce  a  CER  that  can 
be  used  to  estimate  spare  parts  cost. 

The  failure  of  the  data  to  produce  a  usable  CER  is  not 
surprising.  ^he  funding  for  spare  parts  comes  from 
operational  funds.  For  this  reason,  the  cost  of  supplies  is 
a  measure  of  the  availability  of  money  as  well  as  the  cost  of 
replacement  repair  parts. 
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TABLE  4.2 


MULTIPLE  REGRESSION  OF  INDEPENDENT  VARIABLES 
AGAINST  THE  DEPENDENT  VARIABLE  SUPPCOST 


RUN  1 

Model  fitting  results  for:  SUPPCOST 


Independent 

variable 

coefficient 

std.  error  t 

-value 

sig. level 

CONSTANT 

28435.122998 

2.039639E41 

1.3941 

0.1708 

PGMHRS 

-8.033911 

6.666066 

-1.2052 

0.2350 

ACQCOST 

0.002245 

0.001477 

1.5196 

0.1363 

COMCOST 

-0.049499 

0.110375 

-0.4485 

0.6562 

UTILCOST 

0.270332 

0.205072 

1.3182 

0.1947 

R-SQ.  (ADJ.)  =  0.3421  SE=  38015.056957  MAE=  20951.79413 
Previously:  0.0000  0.000000  0.000000 

46  observations  fitted. 

RUN  2 

Model  fitting  results  for:  SUPPCOST 


Independent 

variable  coefficient  std.  error  t-value  sig. level 


CONSTANT 

25588.820008 

1.919855E4 

1.3329 

0.1898 

PGMHRS 

-7.789751 

6.580307 

-1.1838 

0.2432 

ACQCOST 

0.002071 

0.001412 

1.4669 

0.1499 

UTILCOST 

0.232652 

0.185282 

1.2557 

0.2162 

R-SQ.  (ADJ.)  =  0.3546  SE= 
Previously:  0.3421 

46  observations  fitted. 

37651.779438 

38015.056957 

MAE=  20686.445263 
20951.794135 
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C.  DEVELOPMENT  OP  CER  USING  NTSC  DATA 

The  data  provided  by  Naval  Training  Systems  Center  (NTSC) 
consists  of  the  successful  bid  for  each  of  eight  different 
simulator  maintenance  contracts.  These  bids  included  the  cost 
of  five  years  of  maintenance  t-oth  with  and  without  supply 
support  making  it  most  applicable  to  the  questions  at  hand. 

The  results  of  a  simple  regression  of  COMS  cost  without 
supply  support  (COMWOSUP)  as  the  independent  variable  and  cost 
of  supply  support  (SUPPLY)  as  the  dependent  variable  were  very 
good.  The  t-value  of  the  slope  was  9.88510,  compared  to  a 
tabled  value  of  2.998  for  7  degrees  of  freedom  with  an 
alpha=.01,  indicating  the  significance  of  the  coefficient. 
The  R-squared  value  is  a  respectable  94.21.  This  is  probably 
not  good  enough  for  budgeting  purposes,  but  it  is  accurate 
enough  to  be  used  for  the  comparison  of  contracting 
strategies.  The  resulting  CER  is: 

Cost  of  Supply  Support=  -37012.8  +.1565*COMS  Cost  +e  (4.1) 

It  should  be  noted  that  the  cost  for  supply  support  and  the 
COMS  cost  are  for  a  full  five  year  contract. 

The  problem  with  this  CER  is  that  it  is  now  necessary  to 
come  up  with  a  method  to  determine  an  independent  estimate  for 
the  cost  of  a  COMS  contract  before  it  can  be  useful. 
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D.  DEVELOPMENT  OF  CER  FOR  COMS  CONTRACT  COST 

The  regression  analysis  procedures  were  repeated  with 
COMCOST  as  the  dependent  variable-  The  results  were  similar 
to  the  results  of  the  regression  with  SUPPCOST  above.  The  t- 
values  were  good  for  some  of  the  regressions,  but  the  R- 
squared  values  were  not  nearly  high  enough  to  iiidicate  a 
quality  CER.  The  multiple  regression  for  COMCOST  was  not 
successful  either.  Stepwise  regression  also  produced  low  R- 
squared  values,  again  indicating  a  less  thc.n  desirable  CER. 
Hence,  no  CER  exists  in  the  FASO  data  set  for  the  cost  of  a 
COMS  contract. 

E.  SUMMARY  OF  REGRESSION  RESULTS 

Regression  has  yielded  a  CER  for  the  cost  of  supply 
support.  That  CER  uses  the  cost  of  the  COMS  contract  as  an 
independent  variable.  Regression  failed  to  produce  a  usable 
CER  for  the  cost  of  a  COMS  contract:,  so  another  method  must  be 
utilized  to  provide  that  cost.  In  the  next  chapter,  such  an 
alternate  procedure  will  be  demonstrated. 
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V.  SIMULATION  MODEL 


The  failure  of  the  previous  data  sets  to  produce  a  CER 

provides  an  opportunity  to  explore  the  modeling  of  engineering 

standards  to  provide  a  method  of  cost  estimation. 

An  old  proverb  states  that  a  journey  of  a  thousands  miles 
begins  with  but  a  single  step.  To  paraphrase,  a  simulation 
analysis  begins  with  a  single  step  and  must  proceed 
methodically  through  several  steps  if  useful  results  are  to 
be  derived  and  implemented.  But  before  beginning  any 
significant  journey,  it  is  important  to  know  something  about 
the  route  and  noteworthy  milestones  alo’-^  he  way. 


Simulation  analysis  is  a  descriptive  modeling  technique.  As 
such,  it  does  not  provide  the  explicit  problem  formulation 
and  solution  steps  which  are  an  integral  part  of 
optimization  models,  such  as  linear  programming. 
Consequently,  one  must  specify  in  some  detail  a  procedure 
for  the  development  and  use  of  simulation  models  to  assure 
successful  outcomes  from  their  application.  [Ref.  3:pp.  14- 
15] 

The  process  utilized  in  this  thesis  to  insure  accurate 
results  consists  of 

1)  A  definition  of  the  problem  and  understanding  of  the 
process  to  be  simulated  so  a  good  model  formulation  can 
be  constructed; 

2)  A  search  for  data  and  analysis  of  that  data  to  determine 
its  characteristics; 

3)  Development  of  a  model  to  simulate  the  maintenance 
process ; 

4)  Verification  and  validation  of  the  model's  code,  logic 
and  results; 

5)  Utilizing  simulation  results  to  determine  the  cost  of  a 
simulator  maintenance  contract. 
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A.  PROBLEM  DEFINITION  AND  FORMULATION 

Problem  formulation  is,  in  a  sense,  the  most  important 
step  in  a  simulation  analysis.  Appropriate  solutions  to 
inappropriately  formulated  problems  cannot  be  achieved. 
But  before  a  problem  can  be  formulated,  it  must  be 
identified  or  found.  Problem  finding  is,  in  reality, 
choosing  from  among  several  problems  which  are  competing 
for  the  same  resources.  Criteria  for  selection  include 
technical  economic,  and  political  feasibility,  and  the 
perceived  urgency  for  a  solution.  Problem  selection  can 
have  a  significant  impact  on  the  ultimate  success  of  the 
analysis  and  the  implementation  of  the  results.  [Ref.  3:pp 
16-17] 

In  order  to  formulate  this  simulation  model,  it  is  first 
necessary  to  have  an  understanding  of  the  process  to  be 
simulated  and  define  the  problem  to  be  solved. 

1.  The  Simulator  Maintenance  Process 

There  are  two  types  of  simulator  maintenance, 
scheduled  maintenance  and  nonscheduled  maintenance.  Scheduled 
maintenance  is  done  as  required  by  the  preventive  maintenance 
schedule,  is  not  variable,  and  will  therefore  not  be  part  of 
the  model. 

Nonscheduled  maintenance  occurs  as  the  result  of  a 
component  failure  or  system  breakdown.  These  component 
failures  or  system  breakdowns  can  be  discovered  and  reported 
in  one  of  three  scenarios.  First,  maintenance  personnel  can 
find  a  discrepancy  during  scheduled  maintenance  and  report  it. 
Second,  a  pilot  can  find  a  discrepancy  during  the  conduct  of 
a  training  mission  and  not  report  it  until  the  end  of  the 
trainer  period.  Finally,  a  pilot  can  find  a  discrepancy 
during  the  conduct  of  a  training  mission  that  degrades  the 
performance  of  the  simulator  below  that  required  for  the 
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completion  of  the  mission,  and  the  pilot  reports  that 
discrepancy  immediately.  All  three  of  these  scenarios  are 
called  breakdown  events  and  will  be  assumed  to  occur  randomly. 

After  a  discrepancy  has  been  reported,  it  is  assigned 
to  a  technician  for  repair.  This  repair  takes  a  variable 
length  of  time  to  complete.  The  completion  of  the  repair 
makes  the  simulator  available  for  use  and  is  called  a  repair 
event.  The  simulator  is  scheduled  for  use  and  the  process 
starts  over  again. 

2.  Problem  Statement 

To  answer  the  questions  posed  in  this  thesis,  it  is 
necessary  to  determine  how  many  man  hours  are  required  to 
maintain  an  aviation  simulator.  This  will  be  converted  to  the 
number  of  man  years  required  for  maintenance,  and  used  to 
determine  a  cost  estimate. 

B.  SIMULATION  MODEL  DEVELOPMENT 

The  problem  of  determining  the  number  of  man  hours 
required  to  maintain  an  aviation  simulator  can  be  answered  by 
developing  a  Fortran  language  computer  model  to  simulate  the 
maintenance  of  an  aviation  simulator  suite.  This  model  should 
be  easily  modifiable  to  fit  simulator  suites  consisting  of  any 
number  of  elements.  An  element  is  a  segment  of  a  simulator 
suite  designed  for  a  specific  training  mission  (i.e., 
operational  flight  trainer  (OFT) ,  weapons  system  trainer 
(WST) ,  etc) . 
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Such  a  model  was  developed  for  this  thesis.  It  utilizes 
a  main  program  to  initialize  the  parameters  and  perform  the 
simulation  by  calling  a  variety  of  subroutines.  These 
subroutines  are  the  initialization  subroutine,  the  time 
advance  subroutine,  the  breakdown  event  subroutine,  the  repair 
event  subroutine,  the  report  generator  subroutine,  and  random 
variable  generators  for  the  exponential,  uniform,  and  normal 
distributions.  The  failure  and  repair  subroutines  are 
constructed  from  a  block  of  code  with  one  block  for  each 
simulator  element.  This  enables  the  model  to  be  expanded  or 
contracted  to  meet  the  specifications  of  any  simulator  suite. 
At  the  completion  of  the  simulation,  the  report  generator 
prints  out  the  parameters  used  for  the  failure  and  repair 
distributions,  the  number  of  maintenance  personnel,  maximum 
queue  length,  element  operational  ready  rate,  number  of 
maintenance  actions  that  took  longer  than  three  hours  to 
repair,  simulation  run  time,  and  total  number  of  maintenance 
actions  that  occurred  during  the  simulation.  A  complete 
listing  of  the  program  is  included  in  Appendix  A. 

1.  How  to  Use  the  Model 

a.  The  first  step  is  to  load  the  program  into  a  computer 
with  a  Fortran  compiler. 

b.  Modify  the  parameters  for  the  failure  and  repair 
distributions  to  reflect  those  of  the  simulator  suite 
being  simulated.  An  analysis  of  these  distributions  for 
the  example  considered  in  this  thesis  is  in  the  next 
section. 
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c.  Modify  the  number  of  elements  in  the  program  to  reflect 
the  physical  composition  of  the  simulator  suite  of 
interest.  For  the  purposes  of  this  thesis,  the  eight 
elements.  Operational  Flight  Trainer  (OFT)  1  and  2, 
Weapons  System  Trainer  (WST)  1  and  2,  Positional  Type 
Trainer  Tactics  (PTT)  1  and  2,  and  Positional  Type 
Trainer  Cockpit  (PTC)  1  and  2  have  been  used. 

d.  Estimate  the  number  of  maintenance  personnel  required  to 
maintain  the  simulator  suite,  and  enter  that  number  for 
the  value  of  NMTPER.  Then,  determine  the  length  of  time 
the  simulation  is  to  cover  and  enter  that  value  for 
FEL(17) . 

e.  Compile  and  run  the  model. 

f.  Note  the  results  of  interest,  in  this  case  operational 
ready  rate  and  maximum  repair  cpjeue  length,  and  modify 
the  number  of  maintenance  personnel  up  or  down  until  the 
desired  results  are  obtained. 

g.  Determine  the  number  of  maintenance  personnel  required 
to  perform  the  contract.  This  is  done  by  multiplying 
the  number  of  maintenance  personnel  needed  to  meet 
mission  requirements  by  3  to  cover  a  3  shift  24  hour 
day,  and  adding  the  appropriate  number  of  overhead 
personnel.  Overhead  personnel  include,  but  are  not 
limited  to,  site  supervisor,  supply  clerk,  maintenance 
administration  clerk,  personnel  to  perform  the  required 
preventive  maintenance  (PM) ,  and  safety  observers  as 
required  by  Navy  safety  regulations. 

h.  Multiply  the  number  of  man-years  required  to  maintain 
the  simulator  by  the  NTSC-provided  cost  for  acquiring 
one  man-year  on  contract,  $45,000,  to  get  a  cost 
estimate  for  the  COMS  contract. 


C.  DATA  SEARCH  AND  ANALYSIS 

Data  used  in  determining  the  parameters  to  be  utilized  by 
this  model  were  obtained  from  an  actual  P-3  simulator  site. 
Failure  rate  data  were  obtained  for  a  227  day  period  of  time, 
and  time  to  repair  data  were  obtained  for  a  different  90  day 
period  of  time.  The  different  periods  of  time  for  the  two 
data  sets  were  driven  by  the  way  the  records  were  kept  by  the 
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maintenance  contractor.  These  data  were  analyzed  to  determine 
if  the  failure  rate  and  repair  times  demonstrated  the  behavior 
of  a  known  probability  distribution. 

1.  Analysis  of  Data 

a.  Interarrival  Time 

The  periodic  occurrence  of  an  event  such  as  the 
breakdown  of  an  aviation  simulator  is  best  described  in  terms 
of  the  time  between  breakdown  events.  This  period  of  time  is 
called  the  interarrival  time.  The  data  used  to  determine  the 
interarrival  times  were  taken  from  the  maintenance  action 
forms  (MAF)  that  were  filled  out  by  the  person  discovering  the 
malfunction  and  submitted  to  maintenance  personnel.  Once  the 
failure  times  were  entered  into  a  data  bank,  they  were 
processed  by  a  simple  Fortran  program  to  convert  them  into 
interarrival  times. 

The  analysis  of  the  resultant  interarrival  times 
was  done  utilizing  the  IBM  program  GRAFSTAT.  The  "Fitting 
Probability  Distributions”  option  of  the  "Analysis  of 
Univariate  Distributions"  section  of  the  Grafstat  menu  was 
used.  This  option  utilized  the  fitted  histogram/density  and 
probability  plots  to  analyze  the  data  and  determine  its 
parameters.  It  was  initially  anticipated  that  the 
interarrival  times  would  be  exponentially  distributed. 

The  data  used  for  this  thesis  turned  out  to  be 
very  ill  behaved.  The  data  did,  in  fact,  fail  to  conform  to 
any  known  distribution,  and  it  defied  initial  attempts  at 
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analysis.  It  was  therefore  necessary  to  attempt  data 
transformations  in  order  to  get  any  usable  information. 

(1)  Transformation.  Transformation  is  the 
systematic  modification  of  data  to  make  it  conform  to  a 
distribution.  The  most  utilized  of  these  modification 
procedures  are  the  logarithmic  and  exponential 
transformations.  The  process  is  a  simple  one;  the  data  is 
raised  to  a  certain  power,  or  its  logarithm  is  determined,  and 
the  resulting  data  set  is  analyzed.  The  exponential 
transformation  proved  to  be  useful  and  is  represented  by 
equation  5.1, 

Y®  =  X.  (5.1) 
Y  is  the  data,  a  is  the  power  to  which  the  data  is  raised, 
and  X  is  a  random  variable  that  can  be  fitted  to  some 
distribution. 

Unfortunately,  even  after  transformation,  much 
of  the  data  continued  to  defy  analysis.  Probability  plots  of 
the  raw  data  showed  that  much  of  it  was  bimodal  or  lacked  good 
discipline.  Neither  of  these  problems  came  as  a  total 
surprise  because  of  the  inherent  characteristics  of  the 
failure  process. 

(2)  Problems.  A  bimodal  situation  is  created  when 
data  points  on  a  data  set  come  from  two  or  more  distributions. 
That  is  to  say,  in  this  case,  that  the  failure  rate  of  the 
main  computer  is  different  than  that  of  the  motion  system, 
which  is  different  from  that  of  the  cockpit  indications,  etc. 
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When  all  the  types  of  failures  are  thrown  into  one  set  of 
numbers,  it  can  create  many  problems  for  the  analyst. 

The  data  problems  appear  to  be  the  result  of 
the  process  by  which  the  data  were  recorded.  Many  times  a 
discrepancy  will  be  discovered  during  a  training  evolution, 
but,  since  it  does  not  degrade  training,  it  is  not  recorded 
until  the  end  of  the  exercise.  Sometimes  several 
discrepancies  will  be  identified  at  essentially  the  same  time 
(the  completion  of  the  trainer  time  period) .  There  is  also  a 
high  probability  that  several  malfunctions  were  discovered  by 
maintenance  personnel  and  corrected  without  any  paper  work  to 
record  the  malfunction  or  the  repair.  All  of  these  types  of 
data  excursions  cause  the  data  to  be  ill  behaved  and  increase 
the  difficulty  of  analysis. 

(3)  Solution.  To  counter  these  problems,  several 
algorithms  for  the  performance  of  the  required  analysis  were 
devised  and  tested.  The  one  that  produced  the  best  results  is 
described  as  follows. 

The  first  step  of  this  algorithm  was  to 
uniformly  distribute  all  blocks  of  discrepancies  obviously 
written  at  the  end  of  a  training  evolution.  This 
redistribution  was  made  over  the  last  half  of  a  three  hour 
training  period.  The  resulting  data  were  then  analyzed  as 
described  above  in  two  different  ways.  First,  each  element 
(OFTl,  0FT2 ,  etc.)  was  individually  analyzed.  Second,  the 
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elements  were  combined  by  type  (OFT,  WST,  etc.)  and  analyzed 
again. 

The  results  of  this  analysis  are  shown  in 
TABLE  5-1.  For  each  element  and  type,  the  second  column  shows 
exponent  that  was  used  to  transform  the  interarrival  times. 
The  third  column  shows  the  value  of  the  parameter  of  an 
exponential  distribution  which  was  fit  to  the  transformed 
interarrival  times  for  that  element  and  type.  The  fourth 
column  shows  the  significance  level,  using  a  Kolmogorov- 
Smirnov  statistic,  of  the  test  of  the  hypothesis  that  the 
empirical  distribution  of  the  transformed  interarrival  times 
fit  a  theoretical  exponential  distribution  with  the  parameter 
shown  in  the  column  three. 


TABLE  5.1  PARAMETRIC  ANALYSIS  OF  THE  INTERARRIVAL  TIMES  FOR 
SIMULATOR  ELEMENTS 


K.-S. 


ELEMENT 

TRANSFORMATION 

EXPONENTIAL 

SIGNIFICANCE 

EXPONENT 

PARAMETER 

LEVEL 

OFTl 

.84 

51.59 

.9198 

0FT2 

.40 

4.39 

.3131 

OFT 

.  88 

18.52 

.7591 

WSTl 

.  64 

15.77 

.8061 

WST2 

.46 

5.67 

.9039 

WST 

.70 

22.64 

.9057 

PTTl 

.  54 

8.86 

.5918 

PTT2 

.  62 

11.15 

.  6599 

PTT 

.77 

28.05 

.9619 

PTCl 

.  69 

33.65 

.8591 

PTC  2 

.  36 

3.06 

.  1296 

PTC 

.73 

21.99 

.9256 

2^ 


The  most  significant  results  were  then  used  in 
accordance  with  the  following  rules: 

1)  If  analysis  resulted  in  good,  significant 
parameters  for  both  elements  of  a  type  of 
trainer,  those  individual  parameters  were  used  in 
the  simulation. 

2)  If  analysis  resulted  in  very  significant  results 
coming  from  only  one  element  of  a  type,  then  that 
parameter  was  used  for  both  elements  in  the 
simulation. 

3)  If  the  most  significant  results  came  from  the 
combined  elements  of  type  data,  the  parameter  for 
the  combined  data  was  used,  as  is,  for  both 
elements  in  the  simulation. 


In  the  analysis  of  the  failure  data,  rule  2 
was  applied  to  OFTl  and  0FT2  (OFTl  data  being  used) .  Rule  3 
was  applied  to  WSTl,  WST2 ,  PTTl,  PTT2 ,  PTCl,  and  PTC2 .  The 
parameters  used  for  interarrival  times  in  the  simulation  are 
presented  in  TABLE  5.2.  The  effects  of  these  choices  of  the 
parameters  on  the  overall  model  will  be  shown  in  Section  D 
below. 


TABLE  5.2  INTERARRIVAL  TIME  PARAMETERS  USED  IN  THE 
SIMULATION  MODEL 


ELEMENT  TRANSFORMATION 

EXPONENT 


EXPONENTIAL 

PARAMETER 


OFTl 

.84 

51.59 

OFT2 

.  84 

51.59 

WSTl 

.  70 

22.64 

WST2 

.70 

22.64 

PTTl 

.77 

28 . 05 

PTT2 

.  77 

28.05 

PTCl 

.73 

21.99 

PTC2 

.73 

21.99 

30 


b.  Maintenance  Repair  Time 


The  length  of  time  it  takes  a  maintenance 
technician  to  repair  an  aviation  simulator  system  discrepancy 
is  called  the  maintenance  repair  time.  The  maintenance  repair 
times  are  taken  directly  from  the  portion  of  the  MAF  that  is 
filled  out  by  the  technician  after  he  has  finished  the  repair 
process. 


The  analysis  of  the  maintenance  repair  time  data 
was  conducted  identically  to  that  of  the  interarrival  time 
data,  and  the  results.  .  ^re  similar.  The  data  turned  out  to  be 
equally  as  ill  behdved  and  the  transformation  process  was 
again  required  to  obtain  usable  results. 

The  decision  process  was  the  same  as  well,  and 
TABLE  5.3  illustrates  the  transformation  exponents  and  mean 


service  time  parameters  used  in  the  model. 


TABLE  5.3 

MEAN  SERVICE  TIME 
SIMULATION  MODEL 

PARAMETERS 

USED  IN  THE 

ELEMENT 

TRANSFORMATION 

MEAN  SERVICE 

STANDARD 

EXPONENT 

TIME 

DEVIATION 

OFTl 

.  33 

1.14 

.22 

0FT2 

NATURAL  LOG 

2.08 

.  36 

WSTl 

.  32 

1.46 

.43 

WST2 

.  32 

1.46 

.43 

PTTl 

.50 

2.23 

.74 

PTT2 

.  50 

2 .23 

.74 

PTCl 

.  14 

1.05 

.07 

PTC2 

.  14 

1.05 

.  07 
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D.  VERIFICATION  AND  VALIDATION  OF  THE  MODEL 


Model  verification  and  validation  actually  is  concerned  with 
three  models:  a  conceptual  model,  a  logical  model,  and  a 
computer  model.  ...Verification  focuses  on  internal 
consistency  of  a  model,  whole  validation  is  concerned  with 
the  correspondence  between  the  model  and  reality.  [Ref.  3:p. 
27] 

Verification  of  the  model  was  a  continuing  process 
throughout  its  development.  As  they  were  identified,  code  and 
logic  errors  were  rectified  until  the  model  ran 
satisfactorily. 

To  validate  the  simulation  model,  the  parameters  derived 
by  the  analysis  of  available  maintenance  data  were  inserted 
into  the  simulation  model.  The  model  was  set  to  run  for  6000 
hours  which  represents  the  number  of  hours  of  operation 
provided  for  in  one  year  of  operation  (364  days  -  104  weekend 
days  -  10  govt  holidays  *  24  hours/day  =  6000  hours)  .  The 
program  was  run  for  two  maintenance  personnel  (minimum  number 
of  personnel  required  by  safety  regulations  for  maintenance  of 
energized  equipment  with  high  voltage  components)  and 
increased  by  one  person  until  the  maximum  queue  size  was  0. 
The  detailed  results  are  shown  in  Appendix  B. 

It  appears  that  four  maintenance  personnel  are  required  to 
keep  the  simulators  up  without  any  delay  time  due  to  non¬ 
availability  of  maintenance  personnel.  The  preventive 
maintenance  (PM)  manuals  indicate  that  it  takes  3.15  people  to 
do  the  required  PM. 

All  the  information  required  to  determine  the  number 
of  maintenance  personnel  to  maintain  this  simulator  suit  is 
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now  known.  That  number  is  18.15.  TABLE  5.4  illustrates  the 
determination  of  this  number. 

The  simulator  in  question  was  being  maintained  by  18 
people  at  the  time  the  data  was  gathered.  The  model  therefore 
appears  to  be  validated  as  a  good  simulation  of  aviation 
simulator  maintenance. 

TABLE  5.4  NUMBER  OF  PERSONNEL  REQUIRED  TO  MAINTAIN  THE  P-3 
SIMULATOR  MODELED 


BILLET 


NUMBER  REQUIRED 


SIGHT  SUPERVISOR  1 

SUPPLY  CLERK  1 

MAINTENANCE  ADMINISTRATION  CLERK  1 

PM  PERSONNEL  3.15 

MAINTENANCE  PERSONNEL _ 12 


TOTAL 


18.15 


E.  UTILIZATION  OF  MODEL  RESULTS 

The  solution  to  the  problem  of  determining  the  cost  of 
maintenance  and  supply  support  for  aviation  simulators  is  at 
hand.  Simple  computations  are  all  that  is  required  to 
estimate  those  costs. 

The  cost  of  a  contract  requiring  18.15  men  to  maintain 
simulators  is  computed  by  multiplying  18.15  by  $45,000  (the 
approved  estimate  for  the  cost  of  one  man  year  of  an  average 
simulator  technician).  That  cost  is  $816,750  for  one  year  of 
maintenance  service.  This  is  within  5%  of  the  cost  of  one 
year  of  service  as  delivered  by  the  present  contractor. 
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Now  the  developed  COMS  cost  can  be  used  in  the  developed 
CER  (equation  4.1)  to  determine  the  cost  for  contractor  supply 
support  for  this  simulator  suite.  This  cost  is  estimated  to 
be  $127,089  per  year.  This  also  compares  favorably  with  the 
cost  of  the  contract  now  in  effect. 

Therefore,  one  method  of  determining  the  cost  of  supply 
support  provided  by  a  COMS  contractor  has  been  determined. 
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VI.  AREAS  FOR  FURTHER  RESEARCH 


Presently  the  contractor  supplied/Navy  supplied  mix  for 
spare  parts  is  determined  by  what  is  in  the  simulator  supply 
room  at  the  time  the  contractor  takes  over  maintenance  of  the 
simulator.  An  inventory  is  conducted,  and  the  contractor  is 
responsible  for  maintenance  of  that  inventory  until  the 
completion  of  the  contract  period.  While  this  is  a  step  in 
the  right  direction,  there  are  ways  of  determining  a  more  cost 
efficient  spare  parts  mix  that  considers  all  non  trainer 
peculiar  parts. 

The  best  way  to  determine  the  most  cost  efficient  parts 
mix  is  to  construct  a  linear  optimization  program  that 
minimizes  the  cost  of  parts  support.  A  proposed  basic 
formulation  for  th  .t  optimization  program  is: 

INDEX 

i  =  PART  NUMBER 
PARAMETERS 

CNA  =  COST  FOR  NAVY  ACQUISITION  OF  PART 

CNS  =  COST  FOR  NAVY  TO  STORE  PART 

CNH  =  COST  FOR  NAVY  TO  HANDLE  PART 

CCA  =  COST  FOR  CONTRACTOR  ACQUISITION  OF  PART 

CCH  =  COST  FOR  CONTRACTOR  TO  HANDLE  PART 

NR  =  NUMBER  OF  PARTS  REQUIRED  FOR  LIFE  CYCLE  OF 
SIMULATOR 
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VARIABLES 


Xj  =  NUMBER  OF  PARTS  NAVY  TO  BUY 
Yj  =  NUMBER  OF  PARTS  j  CONTRACTOR  TO  BUY 
FORMULATION 

MINIMIZE 

Ej  CNAj*Xj  +  CNSj*Xj  +  CNHj*Xj  +  CCA^*Y^  +  CCHj*Y,- 
SUBJECT  TO 

Xj  +  Yj  =  NRj 

This  basic  formulation  is  simplistic  by  design.  It  has  no 
constraints,  other  than  the  parts  constraint,  to  allow  for 
tailoring  to  any  specific  situation.  For  example,  if  there  is 
a  limited  amount  of  operating  funds  (O&MN)  the  following 
constraint  can  be  added: 

Ej  CNAi*Xi  <  FUNDS  AVAILABLE 
Likewise,  if  there  is  a  paucity  of  funds  available  for 
contractor  supply  support,  the  following  constraint  can  be 
added: 

Ej  CCAj*Yj  +  CCHj*Yj  <  FUNDS  AVAILABLE 
Values  for  Xj  can  be  assigned  for  those  parts  that  are 
trainer  peculiar  and  have  to  be  purchased  by  the  Navy  at  the 
time  of  acquisition  of  the  simulator. 

The  ability  to  tailor  this  formulation  to  any  simulator 
suite  is  limited  only  by  reality.  If  a  constraint  can  be 
identified  and  linearly  represented,  it  can  be  included  in  the 
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formulation.  The  only  real  limitations  are  the  availability 
of  data  and  the  format  of  that  data. 

Data  were  not  available  to  support  this  model  at  the  time 
the  research  was  done  for  this  thesis.  The  hard  copy  data 
that  was  acquired  for  a  small  subsystem  did  not  have  complete 
cost  data  for  over-the-counter  (contractor)  purchase  of  most 
of  the  parts.  A  recommended  format  for  the  data  to  be 
optimized  by  this  formulation  is: 

PART  COST 

INDEX  ACQUISITION  HANDLING  STORE 

NUMBER  NSN _ NAVY  CNTR  NAVY  CNTR  NAVY 

1  1234567890  145.00  128.00  1.50  12.50  13.40 

Once  the  data  have  been  put  into  the  right  format,  the 
formulation  and  its  associated  data  can  be  an  input  into  any 
matrix  generator/linear  solver.  The  resultant  solution  of  the 
linear  set  of  equations  will  provide  the  optimal  mix  of 
contractor/Navy  acquisition  of  spare  parts 

Unfortunately,  the  absence  of  usable  data  prevents  the 
verification  of  this  optimization  model.  However,  the 
formulation's  simplicity,  coupled  with  its  sound  theoretical 
foundation,  would  lend  an  air  of  credibility  to  the  results 
obtained. 
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VII.  SUMMARY  AMD  RECOMMENDATIONS 


A.  SUMMARY 

The  results  of  the  analysis  in  this  thesis  show  that  the 
aircraft  simulator  cost  elements  studied  do  not  correlate  well 
enough  to  allow  for  successful  linear  regression  analysis. 
This  is  most  likely  caused  by  pressure  applied  to  the  level  of 
spending  on  repair  parts  from  sources  other  than  the 
maintenance  process.  Examples  are  requirements  to  expend  all 
funding  during  the  quarter  received  and  the  necessity  to 
transfer  funds  between  commands  to  meet  priority  requirements 
in  other  areas. 

However,  it  has  been  shown  that  there  is  a  predictable 
relationship  between  the  cost  of  a  maintenance  contract  and 
the  cost  for  providing  supply  support  for  that  contract.  This 
relationship  is  derived  from  a  relatively  small  set  of  data 
drawn  from  eight  successful  bidders  for  simulator  maintenance 
contracts.  The  accuracy  of  this  model  is  excellent,  with  a 
high  R-squared  value  (94.21).  Utilization  of  this  model  will 
provide  reasonably  accurate  government  estimates  for 
comparison  with  the  bids  from  commercial  contractors. 

A  maintenance  process  simulation  model  was  created  to 
provide  an  independent  estimate  for  the  cost  of  maintaining  a 
simulator  for  one  year.  Its  accuracy  was  validated  by 
analyzing  the  performance  of  an  actual  simulator  and  utilizing 
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those  performance  parameters  in  the  simulation  model.  The 
results  were  impressive.  In  addition  to  its  accuracy,  the 
model,  because  of  its  modular  construction,  is  easily 
adaptable  to  any  specific  simulator  suite. 

This  thesis  has  also  provided  a  formulation  for  the 
determination  of  the  optimal  mix  of  contractor-supplied  versus 
Navy-supplied  repair  parts.  Even  though  it  has  not  been 
verified  or  validated,  this  optimization  formulation  has  great 
potential  for  savings  if  the  time  is  taken  to  put  the  cost 
data  in  the  format  required. 

B .  RECOMMENDATIONS 

The  results  of  this  study  suggest  that  regression  analysis 
and  computer  modeling  can  be  of  significant  value  in  the 
estimation  of  supply  support  costs  and  maintenance  costs  for 
aircraft  simulators.  Continued  utilization  of  these  analysis 
techniques  to  update  the  CERs  will  further  refine  the  accuracy 
of  the  results  obtained. 

The  potential  cost  savings  that  can  be  gained  by 
utilization  of  the  optimization  formulation  developed  herein 
cannot  be  overemphasized.  The  fervent  pursuit  of  a  usable 
data  base  for  all  simulators  in  the  inventory  is  highly 
recommended  and  will  pay  dividends  in  the  future. 
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APPENDIX  A 


TH8MOD 


TNSMOO 

A  simulation  model  written  to  simulate  the  maintenance  requirements  for  an  aviation  simulator  suit 
and  print  out  the  results. 

WRITTEN  BY 

COR  EDUARD  H.  BIGGERS 
TO  SATISFY  THE  REQUIREMENTS  FOR 
A  MASTER  OF  SCIENCE  DEGREE  IN 
OPERATIONS  RESEARCH 
SEPTEMBER  15,  1989 


DEFINITION 

OF 

VARIABLES 

BRKON 

Break  down  subroutine 

CHKIN 

Time  break  down  occurred 

CHKINO 

Time  break  down  added  to  the  queue 

CIOCJC 

Time  keeper  for  simulation 

ELAUMT 

Simulator  element  waiting  for  maintenance 

EXPON 

Exponential  random  variate  generator 

F 

Number  of  bread  downs  lasting  longer  that  3  hours 

FEL 

Future  events  list 

lAT 

Interarrival  time 

IMEVT 

Imminent  event 

INITLZ 

Initialize  subroutine 

J,K,L 

Do- loop  index 

MIT0F1 

Mean  interarrival  time  0FT1 

MITOfZ 

Mean  interarrival  time  OFT2 

MItUSi 

Mean  interarrival  time  UST1 

MITWS2 

Mean  interarrival  time  UST2 

MITPT1 

Mean  interarrival  time  PTT1 

HITPT2 

Mean  interarrival  time  PTT2 

MITPC1 

Mean  interarrival  time  PTC1 

M1TPC2 

Mean  interarrival  time  PTC2 

MO 

Maximum  length  of  queue 

MST0F1 

Mean  service  time  0FT1 

MST0F2 

Mean  service  time  0FT2 

MSTUS1 

Mean  service  time  WST1 

MSTUS2 

Mean  service  time  WST2 

MSTPT1 

Mean  service  time  PTTI 

MSTPT2 

*«* 

Mean  service  time  PTT2 

MSTPC1 

*•* 

Mean  service  time  PTCl 

MSTPC2 

- 

Mean  service  time  PTC2 

NAUTMT 

- 

NLfnber  of  simulator  elements  waiting  maintenance 

NMTBUS 

- 

Nifnber  of  maintenance  personnel  busy 

NMTPER 

- 

Number  of  maintenance  personnel 

NORML 

- 

Normal  random  variate  generator 

NUMEVS 

Nuiber  of  events 

0PRDY1 

- 

Operational  ready  rate  0FT1 

OPROY2 

- 

Operational  ready  rate  0FT2 

OPROY3 

- 

Operational  ready  rate  WST1 

OPROY4 

- 

O^rational  ready  rate  UST2 
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0PRDY5 

0PRDY6 

0PRDY7 

OPROY8 

REPAIR 

RPTGEN 

SGMOn 

SGMOF2 

SGMWS! 

SGMWS2 

SGMPT1 

SGMPT2 

SGMPCI 

SGNPC2 

SVT 

TIMADV 

TNREP 

TREPTM 

TTHINO 


Operational  ready  rate  PTTI 
Operational  ready  rate  PTT2 
Operational  ready  rate  PTC1 
O^rational  ready  rate  PTC2 
Repair  subroutine 
Report  generator  subroutine 
Standard  deviation  OFT1  service  time 
Standard  deviation  OFT2  service  time 
Standard  deviation  WST1  service  time 
Standard  deviation  WST2  service  time 
Standard  deviation  PTTI  service  time 
Standard  deviation  PTT2  service  time 
Standard  deviation  PTC1  service  time 
Standard  deviation  PTC2  service  time 
Service  time 
Time  advartce  subroutine 
Total  nijii>er  of  repair  events 
Total  repair  time 
Total  time  in  queue 


To  utilize  this  model,  first  determine  the  parameters  for  each  of  the  elements  of  the 
simulator.  These  parameters  include  interarrival  time,  mean  service  time  for  repairs,  and  the  standard 
deviation  for  those  repair  times.  Insert  these  values  into  the  main  program.  Modify  the  break¬ 
down/repair  subroutines  to  emulate  the  number  of  elements  in  the  simulator  suit  being  modeled. 
Determine  the  nunber  of  events  and  enter  that  nunber  into  the  main  program. 


Estimate  the  nunber  of  maintenar>ce  personnel  required  to  maintain  the  simulator  at  the  desired 
operational  ready  rate  and  enter  that  value  in  the  main  program. 


Run  the  simulation  program  and  note  the  results. 


Modify  the  nonber  of  maintenance  personnel  until  the  operational  ready  rate  is  within  an 
acceptable  range,  until  the  maximum  queue  size  is  0,  or  until  any  other  required  condition  is  met. 

**««'******•«••••*•******«**•••**•*«•••«**•••*•*•*«*•**♦♦*•••*•**•*«•*«*•«***«*********************** 
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PROGRAM  TNSMOO 


MAIN  PROGRAM 

1}  INITIALIZE  MODEL 

2)  CALLS  TIME-ADVANCE  AND  EVENT  ROUTINES 

3)  CALLS  REPORT  GENERATOR  TO  TERMINATE  SIMULATION  NORMALLY 


INTEGER  F,ELAUMT(100),TNREP 

REAL  MIT0F1,MIT0F2.MITWS1.MITWS2,MITPT1,MITPT2.MITPC1,MITPC2, 

1  MST0F1,MST0F2,MSTUS1.MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 

2  FEL{17),CHKIN{8),CHKINO(100),NORML.TREPTM(8),TTMINQ(8) 

COMMON  /SIM/  MIT0F1,MIT0F2.MITWS1.MITWS2,NITPT1.MITPT2,MITPCl, 

1  MITPC2,MSTOF1 ,MSTOF2,MSTWSl ,MSTWS2.MSTPT1 ,MSTPT2.MSTPC1 .MSTPC2, 

2  SGMOF 1 , SGM0F2. SGMWS1 , SGNWS2, SGMPT t .SGNPT2, SGMPCI , SGNPC2, NAWTMT, 

3  NMTBUS . CHK I N , CHK I NO , TREPTM , MQ , F , TNREP.NMTPER . ELAUMT , TTM I NO 
COMMON  /TIMEKP/  CLOCK. IMEVT.NUMEVS.FEL 

NUMEVS  =  17 

***»*******»******************«**o*****************************»*********i 

ASSIGN  VALUES  TO  INPUT  PARAMETERS 
(THESE  VALUES  COULD  BE  STORED  IN  A  FILE  AND  READ  INTO  THE  PROGRAM) 

***«**«****»«************•***************•*****************•*************' 


MITOF1 

S 

25.799 

MITOF2 

S 

25.799 

MITWS1 

s 

22.64 

MITWS2 

s 

22.64 

MITPT1 

= 

28.053 

MITPT2 

28.053 

MITPCl 

= 

21.99 

MITPC2 

3 

21.99 

MSTOF1 

3 

1.14 

MSTOF2 

3 

2.0841 

MSTWS1 

3 

1.4604 

MSTUS2 

8 

1.4604 

MSTPT1 

3 

2.2272 

MSTPT2 

3 

2.2272 

MSTPCI 

3 

1.051 

MSTPC2 

3 

1.051 

SGMOF 1 

3 

0.2183 

SGMOF2 

3 

0.35516 

SGMUSl 

3 

0.43201 

SCMUS2 

3 

0.43201 

SGMPT 1 

3 

0.74425 

SGMPT2 

3 

0.74425 

SGMPCI 

3 

0.073798 

SGMPC2 

3 

0.073798 

NMTPER 

3 

2 

#««************************************************************i 

CALL  INITIALIZATION  ROUTINE 

*«*****«««*«********************•*****«**********•*****•*******' 

CALL  INITLZ 

*********«**«*«************************************************' 

CALL  TIME-ADVANCE  ROUTINE  TO  DETERMINE  IMMINENT  EVENT  AND 
ADVANCE  CLOCK  TO  THE  IMMINENT  EVENT  TIME. 

30  CALL  TIMADV 

THE  VARIABLE  'IMEVT'  INDICATES  THE  IMMINENT  EVENT 


IMEVT 

3 

1 

IS 

AN 

OFT1 

BREAK -DOWN 

IMEVT 

3 

2 

IS 

AN 

OFT2 

BREAK-DOWN 

IMEVT 

3 

3 

IS 

A 

WST1 

BREAK-DOWN 

IMEVT 

3 

4 

IS 

A 

WST2 

BREAK-DOWN 

IMEVT 

3 

5 

IS 

A 

PTT1 

BREAK-DOWN 

IMEVT 

3 

6 

IS 

A 

PTT2 

BREAK-DOWN 

IMEVT 

3 

7 

IS 

A 

PTCI 

BREAK-DOWN 

IMEVT 

3 

8 

IS 

A 

PTC2 

BREAK-DOWN 

IMEVT 

3 

9 

IS 

AN 

OFT1 

REPAIR 

IMEVT 

3 

10 

IS 

AN 

OFT2 

REPAIR 

IMEVT 

3 

11 

IS 

A 

WST1 

REPAIR 
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IM€VT  =  12  IS  A  UST2  REPAIR 

IMEVT  =  13  IS  A  PTT1  REPAIR 

IMEVT  -  U  IS  A  PTT2  REPAIR 

IMEVT  =  15  IS  A  PTC1  REPAIR 

IMEVT  =  16  IS  A  PTC2  REPAIR 

IMEVT  =  17  IS  THE  END  OF  THE  SIMULATION 

GO  TO  (40, 40, 40, 40, 40, 40, 40, 40, 50. 50. 50, 50, 50, 50, 50, 50, 60), IMEVT 

CALL  THE  APPROPRIATE  EVENT  SUBROUTINE 

40  CALL  BRKDN 
GO  TO  30 

50  CALL  REPAIR 

GO  TO  TIME  ADVANCE  ROUTINE  TO  SEE  IF  SIMULATION  IS  COMPLETE 
GO  TO  30 


WHEN  SIMULATION  IS  OVER  CALL  REPORT  GENERATOR 

i***«*****«************«****************************i 

60  CALL  RPTGcN 
STOP 
END 


SUBROUTINE  INITLZ 

**«**************************' 

INITIALIZATION  ROUTINE 


INTEGER  F,ELAUMT(100),TNREP 

REAL  M1T0F1,MIT0F2,MITUS1,MITUS2,MITPT1,MITPT2,MITPC1.MITPC2, 

1  MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2.MSTPC1,MSTPC2, 

2  FEL(17),CHKIN(8),CHKINO(100),NORML,TREPTM(8),TTMINO(S) 

COMMON  /SIM/  MITOF1,MITOF2,MITUS1,MITWS2,MITPT1,MITPT2,MITPCl, 

1  MITPC2,MSTOF1,MSTOF2,MSTWSl,MSTUS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 

2  SGMOF 1 , SGMOF  2 , SGMWS 1 , SGMUS2 , SGMPT 1 , SGMPT2 . SGMPC 1 . SGMPC2 , N AWTMT , 

3  NMTBUS,CHKIN.CHKINO,TREPTM,Ma,F,THREP,NMTPER.ELAWMT,TTMINQ, 
COMMON/TIMEKP/  CLOCK, IMEVT, NUMEVS.FEL 

’*****•*«***«******«*•*******•*•**•**************•********************** 

INITIALIZE  SIMULATION 

1)  SET  SIMULATION  CLOCK  TO  0 

2)  ASSUME  SYSTEM  IS  EMPTY  AND  IDLE  AT  TIME  0 

3)  INITIALIZE  CUMULATIVE  STATISTICS  TO  0 

*********************************************************************** 

DO  10  I  =  1,8 
TTMINQ(I)  =  0 
TREPTM(I)  =  0 
10  CONTINUE 

CLOCK  =  0.0 
IMEVT  =  0 
NAWTMT  =  0 
NMTBUS  =  0 
TNREP  =  0 
F  =  0 
MO  =  0 

«********«*******«*****************************************************! 

GENERATE  TIMES  FOR  THE  FIRST  BREAKDOWNS  AND  SCHEDULE  THEM  IN 
FEL(1)  -  FEL(8),  SET  FEL(9)  -  FEL(16)  TO  INFINITY  TO  INDICATE 
THAT  REPAIR  IS  NOT  POSSIBLE  WHILE  THE  SYSTEM  IS  UP,  SET 
FEL(17)  EQUAL  TO  LENGTH  OF  TIME  SIMULATION  IS  TO  RUN. 

FEL(I)  *  CLOCK  ♦  EXP0N(MIT0F1))**(1/.84) 

FEL(2)  *  CLOCK  ♦  EXPON(MITOF2))**(1/.84) 

FEL(3)  =  CLOCK  ♦  EXPON(MITWS1))*»(1/.7) 

FEL(4)  =  CLOCK  ♦  EXPON(MITWS2))**<1/.7) 

FEL(5)  =  CLOCK  ♦  EXPON(MITPT1 ))**(1/.77) 

FEL(6)  =  CLOCK  ♦  EXPON(MITPT2))**{I/.77) 

FEL(7)  «  CLOCK  ♦  EXPON(MITPC1))**(1/.73) 

FEL(8)  *  CLOCK  +  EXPON(MITPC2))»*(1/.73) 
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FEL(9)  =  1.0E+30 
FEL(IO)  =  1.0E+30 
FEL(II)  =  1.0E+30 
FEL(12)  =  1.0E+30 
FEL(13)  =  1.0E+30 

FEL(14)  =  1.0E+30 

FEL<15)  =  1.0E+30 
FEL(16)  =  1.0E+30 

FEL(17)  =  6000 
RETURN 


SUBROUTINE  TINADV 

TIME  ADVANCE  ROUTINE  FINDS  NEXT  EVENT  ON  FUTURE  EVENT  LIST 
AND  ADVANCES  CLOCK 


INTEGER  F.ELAUMT(100),TNREP 

REAL  M1T0F1,MIT0F2,MITWS1,MITWS2,MITPT1,M1TPT2,MITPC1,MITPC2, 

1  MST0F1 ,MST0F2,MSTUS1 ,MSTWS2,MSTPT1 ,MSTPT2,MSTPC1 .MSTPC2. 

2  FEL(17),CHKIN(S),CHKINO(100),NORML,TREPTM(8).TTHINO(8) 

COMMCW  /SIM/  MIT0F1,MIT0F2,MITWS1,MITWS2,MITPT1,MITPT2,HITPCl, 

1  MITPC2,MST0F1,MST0F2,MSTWSI,MSTWS2.MSTPT1,MSTPT2,MSTPC1.MSTPC2, 

2  SGMOF 1 , SGM0F2 , SGMWS 1 , SGMWS2 , SGMPT 1 . SGMPT2 . SGMPC 1,SGMPC2,NAUTMT, 

3  NMTBUS . CHK I N , CHK I NQ , TREPTM , MO , F , TNREP, NMTPER , ELAUMT , TTM I NO 
COMMON  /TIMEKP/  CLOCK. IMEVT.NUMEVS.FEL 

FMIN  =  1.0E+29 
IMEVT  =  0 

••***•••*****«••**•***•**•*****••**•••**•••*••••*•**•*****••••*•**«••*•••***•***••**•******••****•*** 

SEARCH  THE  FUTURE  EVENTS  LIST  FOR  THE  NEXT  EVENT 

**************************************************i^************************************************** 

DO  30  I^l.NUMEVS 

IF  (FEL(I).GE.FMIN)  GO  TO  30 
FMIN  =  FEL(I) 

IMEVT  »  I 
30  CONTINUE 

IF  < IMEVT. GT.O)  GO  TO  50 

***************************************************************************************************** 

ERROR  CONDITION  -  FUTURE  EVENTS  LIST  IS  EMPTY 

***************************************************************************************************** 

40  FORMAT(2X,  '•***•  FUTURE  EVENTS  LIST  EMPTY******, 

C/, IX, ******  SIMULATION  CANNOT  CONTINUE  ****•<) 

CALL  RPTGEN 
STOP 


ADVANCE  SIMULATION  CLOCK 

NEXT  EVENT  IS  TYPE  *IMEVT'  WHICH  WILL  OCCUR  AT  TIME  FELCIMEVT) 

****•*•«•**•**••«•*••«**•••••*•••*•**•••••**•*••**•*•***»••**••••••*•*• 

50  CLOCK  =  FELCIMEVT) 

RETURN 

END 

SUBROUTINE  BRKDN 

**•••**•**•*•  •*******************************************«**********tt**^ 

BREAK-DOWN  EVENT  ROUTINE  -  SCHEDULES  THE  BREAK-DOWN  EVENT 
MAINTENANCE,  OR  TO  A  QUEUE  AWAITING  MAINTENANCE.  COMPUTE 
MAINTENANCE  TIME  AND  UPDATE  STATISTICS. 


INTEGER  F.ELAWMTC 100), TNREP 

REAL  MIT0F1,MIT0F2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 

1  MSTOF1 ,MSTOF2,MSTWS1 ,MSTWS2,MSTPT1 ,MSTPT2,MSTPC1 ,MSTPC2, 

2  FELC17},CNKIN(8),CHKINO(100),NORML.TREPTM(8),TTMINO(8} 

COMMON  /SIM/  MIT0F1,MIT0F2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCI, 

1  MITPC2,MSTOF1,MSTOF2,MSTWSl,MSTWS2,MaTPT1,MSTPT2.MSTPC1,MSTPC2, 

2  SGMOF 1 , SGMOF2 , SGMWS  1 , SGMWS2, SGMPT 1 , SGMPT2, SGMPC1 , SGMPC2 , NAWTMT , 

3  NMTBUS, CHKIN, CHKINO, TREPTM, MO, F, TNREP, NMTPER. ELAWMT.TTMINO 
COMMON  /TIMEKP/  CLOCK, IMEVT .NUMEVS, FEL 

SVT  =  0 

************************************************************************ 


MAINTENANCE  PERSONNEL  ARE  AVAILABLE.  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 

***************************************************************************************************** 

NHTBUS  X  NHTBUS  *  ^ 

CHKIN(l)  -  CLOCK 

SVT  =  NORML(MSTOF1,SGMOF1) 

FEL(9)  =  CLOCK  +  SVT*»(1/.33) 

FEL(I)  =  1.0E+30 
GO  TO  26 

****************************************************************************************************** 

MAINTENANCE  PERSONNEL  NOT  AVAILABLE.  PUT  REPAIR  IN  THE  Q 
UPDATE  STATISTICS 

•A*************************************************************************************************** 

25  NAUTMT  =  NAUTMT  +  1 

IF  (NAWTMT.GT.MO)  HO  =  NAUTMT 
IF  <NAWTMT.GT.99)  GO  TO  100 
CHKINO(NAUTMT)  =  CLOCK 
ELAUMT( NAUTMT)  -  1 
FEL(I)  >  1.0E+30 

26  RETURN 

***************************************************************************************************** 

DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  BUSY 

30  IF(NMTBUS.GE.NMTPER)  GO  TO  35 

***************************************************************************************************** 


MAINTENANCE  PERSONNEL  ARE  AVAILABLE.  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 

**«**•*♦*♦**«**«•*«*••*****************•*******•**••**********•**** 

NMTBUS  =  NMTBUS  *  1 

CHKIN(2)  =  CLOCK 

SVT  *  NORML(MSTOF2,SGMOF2) 

FEL(IO)  =  CLOCK  ♦  EXP(SVT) 

FEL(2)  =  1.0E+30 
GO  TO  36 


MAINTENANCE  PERSONNEL  NOT  AVAILABLE.  PUT  REPAIR  IN  THE  0 
UPDATE  STATISTICS 

•*•****•••*•**••••«*••«*••*•**••••••**••*••*•••*••••••***•••**••••*••••***•*•••**•••*••••*•*•****•*** 

35  NAUTMT  :  NAUTMT  *  1 

IF  ( NAUTMT. GT. MO)  MO  =  NAUTMT 
IF  < NAUTMT. GT. 99)  GO  TO  100 
CHK I  NO (NAUTMT)  >  CLOCK 
ELAUMT(HAUTMT)  =  2 
FEL(2)  *  1.0E+30 

36  RETURN 

DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  BUSY 
40  IF(NMTBUS.GE.NMTPER)  GO  TO  45 

****ik***********«**  *******  *************************************************************************** 

MAINTENANCE  PERSONNEL  ARE  AVAILABLE.  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 


NMTBUS  =  NHTBUS  +  1 

CNKINO)  =  CLOCK 

SVT  s  N0RML(MSTUS1,SGMUS1) 

FEL(II)  =  CLOCK  +  SVT**(1/.32) 

FEL(3)  =  1.0E+30 

GO  TO  46 


MAINTENANCE  PERSONNEL  NOT  AVAILABLE,  PUT  REPAIR  IN  THE  Q 
UPDATE  STATISTICS 


45  NAWTMT  =  NAUTMT  +  1 

IF  ( NAUTMT. GT. MO)  MQ  =  NAUTMT 
!r  ,NAWIMT.GT.99)  GO  TO  100 
CHK I  NO (NAUTMT)  =  CLOCK 
ELAUMT( NAUTMT)  =  3 
FEL(3)  =  1.0E+30 

46  RETURN 

*************«**«*************************************i 

DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  BUSY 


50  IF(NMTBUS.GE.NMTPER)  GO  TO  55 

***************************************************************************************************** 

MAINTENANCE  PERSONNEL  ARE  AVAILABLE,  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 

***************************************************************************************************** 

NMTBUS  =  NMTBUS  1 

CHKIN(4)  =  CLOCK 

SVT  =  NORML(MSTUS2,SGMUS2) 

FEL(12)  =  CLOCK  ♦  SVT**{l/.32) 

FEL(4)  =  1.0E+30 
GO  TO  56 

***************************************************************************************************** 

MAINTENANCE  PERSONNEL  NOT  AVAILABLE,  PUT  REPAIR  IN  THE  0 
UPDATE  STATISTICS 

***************************************************************************************************** 

55  NAUTMT  *  NAUTMT  ♦  1 

IF  ( NAUTMT. GT. MO)  MQ  =  NAUTMT 
IF  ( NAUTMT. GT. 99)  GO  TO  100 
CHK I NQ( NAUTMT)  s  CLOCK 
ELAUMT (NAUTMT)  =  4 
FEL(4)  =  1.0E+30 

56  RETURN 

***************************************************************************************************** 

DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  U)SY 

***************************************************************************************************** 

60  IF(NMTBUS.GE.NMTPER)  GO  TO  65 

*************************** ************************************************************************** 

MAINTENANCE  PERSONNEL  ARE  AVAILABLE,  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 


NMTBUS  =  NMTBUS  *  1 

CHKIN(5>  =  CLOCK 

SVT  =  NORML(MSTPT1,SGMPT1) 

FEL(13)  *  CLOCK  ♦  SVT**(V.5) 

FEL(5)  =  1.0E+30 
GO  TO  66 

***********  ****************************************************1 

MAINTENANCE  PERSONNEL  NOT  AVAILABLE,  PUT  REPAIR  IN  THE  Q 
UPDATE  STATISTICS 


65  NAUTMT  =  NAUTMT  «  1 

IF  ( NAUTMT. GT. MO)  MO  =  NAUTMT 
IF  (NAUTMT.GT.99)  GO  TO  100 
CHKINO(NAUTHT)  =  CLOCK 
ELAUMT(NAUTMT)  =  5 
FEL(5)  =  1.0E+30 

66  RETURN 
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DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  BUSY 


70  IF(NMTBUS.GE.NMTPER)  GO  TO  75 

MAINTENANCE  PERSONNEL  ARE  AVAILABLE.  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 


NMTBUS  =  NMTBUS  *  1 

CHKIN(6)  =  CLOCK 

SVT  *  NORML(MSTPT2,SGMPT2) 

FEL(14)  =  CLOCK  ♦  SVT**(1/.5) 

FEL{6)  =  1.0E+30 

GO  TO  76 


MAINTENANCE  PERSONNEL  NOT  AVAiLABLE.  PUT  REPAIR  IN  THE  Q 
UPDATE  STATISTICS 


75  NAWTMT  =  NAWTMT  ♦  1 

IF  ( NAWTMT. GT.MQ)  MQ  =  NAWTMT 
IF  ( NAWTMT. GT. 99)  GO  TO  100 
CHKINO(NAWTMT)  =  CLOCK 
ELAUMTC NAWTMT)  =  6 
FEL(6)  =  1.0E+30 

76  RETURN 

****»*«»  >*»****»»****»»»*»*********«*«*»**«*******»»»*i 

DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  BUSY 

***********************  *******************************1 

80  IF(NMTBUS.GE.NMTPER)  GO  TO  85 


MAINTENANCE  PERSONNEL  ARE  AVAILABLE.  UPDATE  SYSTEM  STATE 
COMPUTE  SERVICE  TIME  FOR  REPAIR 

**************  *************************************************1 

NMTBUS  =  NMTBUS  *  1 

CHKIN(7)  *  CLOCK 

SVT  =  MOfiML<MSTPC1.SGMPC1) 

FEL(15)  *  CLOCK  ♦  SVT**(1/.U) 

FEL<7)  =  1.0E+30 
GO  TO  86 

MAINTENANCE  PERSONNEL  NOT  AVAILABLE,  PUT  REPAIR  IN  THE  Q 
UPDATE  STATISTICS 

**************  *************************************************1 

85  NAWTMT  =  NAWTMT  ♦  1 

IF  ( NAWTMT. GT. MO)  HO  =  NAWTMT 
IF  ( NAWTMT. GT. 99)  GO  TO  100 
CHK I NO( NAWTMT)  =  CLOCK 
ELAWMT( NAWTMT)  =  7 
FEL(7)  =  1.0E+30 

86  RETURN 

******  *********************************************************1 

DETERMINE  IF  MAINTENANCE  PERSONNEL  ARE  ALL  BUSY 

90  IF(NMTBUS.GE.NMTPER)  GO  TO  95 


MAINTENANCE  PERSONNEL  ARE  AVAILABLE,  UPDATE  SYSTEM  STATE  AND 
COMPUTE  SERVICE  TIME  FOR  REPAIR 


NMTBUS  =  NMTBUS  *  1 

CHKIN(8)  =  CLOCK 

SVT  =  MOfiML(MSTPC2,SGMPC2) 

FELC16)  *  CLOCK  ♦  SVT**<1/.14) 

FEL(8)  =  1.0E+30 
GO  TO  96 

***************************************************************1 

MAINTENANCE  PERSONNEL  NOT  AVAILABLE,  PUT  REPAIR  IN  THE  0 
UPDATE  STATISTICS 

***************************************************************! 


r*** 
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95  MAUTMT  =  NAUTMT  +  1 

IF  ( NAUTMT. GT.MQ)  HQ  =  NAUTMT 
IF  (NAUTMT. GT. 99)  GO  TO  100 
CHKINO(NAUTMT)  =  CLOCK 
ELAUMT( NAUTMT)  =  8 
FEL(8)  =  1.0E+30 

96  RETURN 

100  UR1TE{6,105) 

105  FORMAT (3X, TO  MANY  AUAITING  MAINTENANCE  ACTIONS 

C1X, INCREASE  NUMBER  OF  MAINTENANCE  PERSONNEL  ***•*■) 

CALL  RPTGEN 

STOP 

END 

SUBROUTINE  REPAIR 

REPAIR  EVENT  ROUTINE:  NOTE  REPAIR  COMPLETE,  TAKE  NLXT 
BREAK-OOUN  EVENT  OUT  OF  THE  QUEUE,  COMPUTE  REPAIR  TIME, 

SCHEDULE  REPAIR  EVENT. 

INTEGER  F,ELAUHT(100),TNREP 

REAL  MITOF1,MITOF2,MITUS1,MITUS2,MITPT1,MITPT2,MITPC1,MITPC2, 

1  MSTOF1,MSTOF2,MSTUS1,MSTUS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 

2  FEL(17),CHKIN(8),CHKINQ(100),NORHL,TREPTH(8),TTMINQ(8) 

COMMON  /SIM/  MITOF1,MITOF2,MITUS1,MITUS2,MITPT1,MITPT2,MITPCl, 

1  MITPC2,MSTOF1,MSTOF2,MSTUSl,MSTUS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 

2  SGMOF 1 , SGMOF2 , SGMUS 1 , SGMUS2 , SGMPT 1 , SGHPT2, SGMPC1 , SGHPC2 , NAUTMT , 

3  NMTBUS,CHKIN,CHKINO,TREPTM,MQ,F,TNREP,NMTPER,ELAUHT,TTMINQ 
COMMON  /TIMEKP/  CLOCK,  IMEVT ,NLIHEVS, FEL 

***************************************************************************************************** 

SCHEDULE  APPROPRIATE  TYPE  OF  REPAIR 
20  -  OFT1  REPAIR 
30  -  OFT2  REPAIR 
40  •  UST1  REPAIR 
50  -  UST2  REPAIR 
60  -  PTT1  REPAIR 
70  ■  PTT2  REPAIR 
80  -  PTC1  REPAIR 
90  -  PTC2  REPAIR 

••M**********************************************************************-*************************** 

GO  TO  (20,30,40,50,60,70,80,90), IHEVT-8 

***************************************************************************************************** 

UPDATE  CUMULATIVE  STATISTICS 

***************************************************************************************************** 

20  TREPTM(I)  =  TREPTH(I)  +  (CLOCK  -  CHKIN(l)) 

TNREP  =  TNREP  ♦  1 

IF((CLOCK  -  CHKIN(1)).GE.3.0)  F  =  F  ♦  1 

*«*********«***♦#*-*******•***************************•*****************«******«********************** 

CHECK  CONDITION  OF  QUEUE 

c**************************************************************************************************** 

IF(NAUTMT).GE.I.O)  GO  TO  25 

QUEUE  EMPTY,  SERVICE  MAN  BECOMES  IDLE,  NEXT  DEPARTURE  TIME 
EQUAL  TO  INFINITY 

NMTBUS  =  NMTBUS  -  1 
FEL(9)  =  1.0E+30 
GO  TO  28 

****♦*»****#***♦♦♦*****♦»♦ ***•#♦»»*»** **«♦»♦**»»***»»»**»»»»***»*»*♦♦»»**»♦»*****»»**♦*******»***»*** 

QUEUE  CONTAINS  AT  LEAST  ONE  MORE  MAINTENANCE  ACTION,  START  UORK  ON 
THE  NEXT  MAINTENANCE  ACTION  IN  THE  Q,  AMO  tWVE-THE  REST  UP  ONE 
POSITION  IN  THE  STACK 

***************************************************************************************************** 

25  GO  TO  (26,36,46,56,66,76,86,96),ELAUMT(t) 

26  CHKIM(I)  =  CHKINQ(I) 

TTMINQ(I)  =  TTMIHQ(I)  ♦  (CLOCK  -  CHKIHQ(l)) 

DO  27  I  =  1, NAUTMT  -  1 
CHKINOd)  =  CHKINO(dl) 

ELAUMT(I)  =  ELAUMT(I+l) 
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27  CONTINUE 


UPDATE  LENGTH  OF  QUEUE 

«******>*****■■«********«***********' 

NAWTMT  =  NAWTMT  -  1 

SVT  =  NOftHKMSTOFT.SGNOFI) 

FEL(9)  =  CLOCK  ♦  SVT**(l/.33) 

RETURN 


SCHEDULE  NEXT  BREAK-DOWN  FOR  OFT1 


28  lAT  =  EXPON(MITOFI) 

FEL(1)  =  CLOCK  +  lAT**(l/.84) 
RETURN 


UPDATE  CUMULATIVE  STATISTICS 


30  TREPTM(2)  =  TREPTM(2)  ♦  (CLOCK  -  CHKIN(2)) 

TNREP  =  TNREP  +  1 

IF({CLOCK  -  CHKIN(2)).GE.3.0)  F  =  F  ♦  1 

»♦♦»*>***»»»  ***♦*»****> ****************************************************************************** 

CHECK  CONDITION  Ct  QUEUE 

***************************************************************************************************** 

IF(NAWTMT.GE.I.O)  GO  TO  35 

***************************************************************************************************** 

QUEUE  EMPTY,  SERVICE  MAN  BECOMES  IDLE,  NEXT  DEPARTURE  TIME 
EQUAL  TO  INFINITY 


NMTBUS  =  NMT8US  -  1  FEL(IO)  =  1.0E+30 
GO  TO  38 

***************************************************************************************************** 


QUEUE  CONTAINS  AT  LEAST  ONE  MORE  MAINTENANCE  ACTION,  START  WORK 
ON  THE  NEXT  MAINTENANCE  ACTION  IN  THE  Q,  AND  MOVE  ONE  THE  REST 
UP  ONE  POSITION  IN  THE  STACK 

***************************************************************************************************** 

35  GO  TO  (26,36,<;6,56,66,76,86,96),ELAWMT(() 

36  CHKIN<2)  =  CHKINQ(I) 

TTMINQ(2)  =  TTHIN0(2)  ♦  (CLOCK  -  CHKINQ(l)) 

DO  37,  I  =  1, NAWTMT  -  1 
CHKINQ(I)  =  CHKINQ(I+l) 

ELAUHT(I)  =  ELAUMT(I*l) 

37  CONTINUE 


UPDATE  LENGTH  OF  QUEUE 

************************************* 

NAWTMT  =  NAWTMT  -  1 

SVT  =  NOfiML(MSTOF2,SGMOF2) 

FEL(IO)  =  CLOCK  ♦  SVT**(l/.33) 

RETURN 


SCHEDULE  NEXT  BREAK-DOWN  FOR 

*************************************************' 

38  lAT  =  EXPON(MITOF2) 

FEL(2)  =  CLOCK  ♦  IAT**(l/.84) 

RETURN 

****************************  *********************^ 

UPDATE  CUMULATIVE  STATISTICS 

*************************************************' 

40  TREPTM(3)  =  TREPTM(3)  ♦  (CLOCK  -  CHKIN(5)) 
TNREP  =  TNREP  ♦  1 

IF((CLOCK  -  CHKIN(3)).GE.3.0)  F  =  F  *  1 

*************************************************' 

CHECK  CONDITION  OF  QUEUE 


IF(NAWTMT).GE.1 .0)  GO  TO  45 


QUEUE  EMPTY,  SERVICE  MAN  BECOMES  IDLE,  NEXT  DEPARTURE  TIME 
EQUAL  TO  INFINITY 
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NMTBUS  =  NHTBUS  -  1 
FELdI)  =  1.0E+30 
GO  TO 


QUEUE  CONTAINS  AT  LEAST  ONE  MORE  MAINTENANCE  ACTION,  START  WORK  ON 
THE  NEXT  MAINTENANCE  ACTION  IN  THE  Q,  AND  MOVE-THE  REST  UP  ONE 
POSITION  IN  THE  STACK 

45  GO  TO  (26,36,46,56.66.76.86,96),ELAUMT(l) 

46  CHKINO)  =  CHKINOd) 

TTMINO(3)  =  TTMINQ(3)  +  (CLOCK  -  CHKINOd)) 

DO  47  I  =  I.NAWTMT  -  1 
CHKINOd)  =  CHKINQ(dl) 

ELAUHTd)  =  ELAUMTd  +  l) 

47  CONTINUE 


UPDATE  LENGTH  OF  QUEUE 


NAUTMT  =  NAUTMT  -  1 

SVT  =  NORHL(MSTWS1,SGMWSt) 

FELdI)  =  CLOCK  ♦  SVT**d/.32) 

RETURN 

***************************************************************************************************** 

SCHEDULE  NEXT  BREAK-DOWN  FOR  WST1 

*-*******************************«******************************************************************** 

48  I AT  =  EXPON<MITUS1) 

F£L<3)  =  CLOCK  +  IAT**d/.7) 

RETURN 

***************************************************************************************************** 

UPDATE  CUMULATIVE  STATISTICS 

**•**»»»***••••»•*»***••»»***•»********»****»»**•»»»»***»»»»»**»»»»*»»»»»»*»»***»**»*****»**»**»»***• 

50  TREPTM(4)  =  TREPTM(4)  +  (CLOCK  -  CHKIN(4)) 

TNREP  =  TNREP  +  1 

IF((CLOCK  •  CHKIN(4)).GE.3.0)  F  =  F  +  1 

**•*•••*••***•«•**********•**•****•***•*•••***■»»****»***»***•*•***•••****•*•••**••••*•••••••*••*••*•• 

CHECK  CONDITION  OF  QUEUE 

****«•**•******••••••••***••••**•***•*••**•••***•*»*•*•••*•••••**••••****•••**•*•******•*•**•***••••• 

IF(NAUTMT.GE.I.O)  GO  TO  55 

•**•••*•»*«*••**»**••*••***•••**•*••••***•*«•***••********•****»***•*•**•••»••****••••••••**••••*•*•• 

QUEUE  EMPTY,  SERVICE  MAN  BECOMES  IDLE,  NEXT  DEPARTURE  TIME 
EQUAL  TO  INFINITY 

***************************************************************************************************** 

NMTBUS  =  NHTBUS  -  1 
FEL(12)  =  1.0E+30 
GO  TO  58 

***************************************************************************************************** 

QUEUE  CONTAINS  AT  LEAST  ONE  MORE  MAINTENANCE  ACTION,  START  WORK 
ON  THE  NEXT  MAINTENANCE  ACTION  IN  THE  0,  AND  MOVE  ONE  THE  REST 
UP  ONE  POSITION  IN  THE  STACK 

***************************************************************************************************** 

55  GO  TO  (26, 36, 46, 56, 66, 76, 86, 96), ELAUHTd) 

56  CHKIN(4)  =  CHKINOd) 

TTMINQ(4)  =  TTMIN0(4)  +  (CLOCK  -  CHKINOd)) 

DO  57,  I  =  1, NAUTMT  -  1 
CHKINOd)  =  CHKINOd  +  l) 

ELAUHTd)  =  ELAUMT(I  +  l) 

57  CONTINUE 


UPDATE  LENGTH  OF  QUEUE 


NAUTMT  =  NAUTMT  -  1 

SVT  =  M0RML(MSTUS2,SGMUS2) 

FEL(12)  =  CLOCK  ♦  SVT»*d/.32) 
RETURN 

SCHEDULE  NEXT  BREAK-DOWN  FOR  WST2 

*********  t******************************^ 
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58  lAT  =  EXf>0N(mTWS2) 

FEL(4)  =  CLOCK  +  IAT**(l/.7) 
RETURN 

UPDATE  CUMULATIVE  STATISTICS 


QUEUE  CONTAINS  AT  LEAST  ONE  MORE  MAINTENANCE  ACTION,  START  WORK 
ON  THE  NEXT  MAINTENANCE  ACTION  IN  THE  0.  AND  MOVE  ONE  THE  REST 
UP  ONE  POSITION  IN  THE  STACK 

****«***1H^********«**********«**««*************************************< 

75  GO  TO  (26,36,46,56,66,76,86,96),ELAUHT(l) 

76  CHKIH(6)  =  CHKINO(I) 

TTMIN0(6)  =  TTHINQ(6)  ♦  (CLOCK  -  CHKINQd)) 


MAUTMT  =  NAWTMT  -  1 

SVT  X  N0RHL(MSTPT2,SGMPT2) 

FEL(H)  =  CLOCK  ♦  SVT**(l/.5) 

RETURN 


80  TREPTM(7)  =  TREPTM{7)  ♦  (CLOCK  -  CHKIN(7)) 

TNREP  =  TNREP  ♦  1 

(F((CLOCK  -  CHKIN(7)).GE.3.0)  F  =  F  ♦  1 
CHECK  CONDITION  OF  QUEUE 

•**•*******•**•*••*•***•*•••*******•*••***•**••*•*****••••*******•••••••*•••••********••••••*******•• 

IF(NAUTMT).GE.I.O)  GO  TO  85 

***************************************************************************************************** 

QUEUE  EMPTY,  SERVICE  MAN  BECOMES  IDLE.  NEXT  DEPARTURE  TIME 
EQUAL  TO  INFINITY 

*********«***************/•*************************************************************************** 

NMTBUS  X  NMTBUS  -  1 
FELdS)  X  1.0E+30 
GO  TO  88 

***•*•**•*•*•••••**•••••**••*••••••*••**••••**••*•*•••***•••••••***•••••*•*•••*••••••*>•*•***••**•**• 

QUEUE  CONTAINS  AT  LEAST  ONE  MORE  MAINTENANCE  ACTION,  START  WORK  ON 
THE  NEXT  MAINTENANCE  ACTION  IN  THE  Q,  AND  MOVE-THE  REST  UP  ONE 
POSITION  IN  THE  STACK 

85  GO  TO  <26,36,46,56,66,76,86,96),ELAUMT(l) 

86  CHKIN(7)  X  CHKINOd) 

TTMINQ(7)  X  TTMINQ(7)  ♦  (CLOCK  -  CHKINO(l)) 

DO  87  I  X  1, NAWTMT  -  1 
CHKINQ(I)  X  CHKINQ(I-^l) 

ELAWMT(I)  X  ELAWMT(I+l) 

87  CONTINUE 

****************x****************************************************************x**x**************** 

UPDATE  LENGTH  OF  QUEUE 

x**************************************************************************************************** 
NAWTMT  X  NAWTMT  -  1 
SVT  X  NORML(MSTPC1,SGMPC1) 

FEL(15)  X  CLOCK  ♦  SVT*»(l/.U) 

RETURN 

SCHEDULE  NEXT  BREAK-DOWN  FOR  PTT1 


88  lAT  X  EXPON(MITPCI) 

FEL(7)  X  CLOCK  +  IAT**(l/.73) 
RETURN 


UPDATE  CUMULATIVE  STATISTICS 


90  TREPTM(8)  x  TREPTM(8)  ♦  (CLOCK  -  CHKIN(8)) 
TNREP  X  TNREP  +  1 

IF((CLOCK  -  CHKIN(8)).GE.3.0)  F  x  f  ♦  1 
CHECK  CONDITION  OF  QUEUE 


IF(NAWTMT.GE.I.O)  GO  TO  95 
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QUEUE  EMPTY,  SERVICE  NAN  BECOMES  IDLE,  NEXT  DEPARTURE  TIME 
EQUAL  TO  INFINITY 


WRITE  OUTPUT 


WRITE(6,10) 

10  F0RMAT(6X,  'SIMULATION  MODEL  OF  AIRCRAFT  SIMULATOR  MAINTENANCE', 

1///) 

WRITE{6,11)  MIT0F1,MIT0F2,MITWS1,MITUS2,MITPT1,MITPT2,MITPC1, 
IMITPC2,MST0FI,MST0F2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2 


11  FORMATL IX, 'LAMBDA  USED  FOR  OFT1  FAILURE  TIMES  ',F10.2 

1  /IX, 'LAMBDA  USED  FOR  OFT2  FAILURE  TIMES  ',F10.2 

2  /IX, 'LAMBDA  USED  FOR  WSTI  FAILURE  TIMES  ',F10.2 

3  /IX, 'LAMBDA  USED  FOR  WST2  FAILURE  TIMES  ',F10.2 
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1 


12 


13 


4 

/IX, 'LAMBDA  USED  FOR  PTT1  FAILURE  TIMES 

',F10.2 

5 

/IX, 'LAMBDA  USED  FOR  PTT2  FAILURE  TIMES 

',F10.2 

6 

/IX, 'LAMBDA  USED  FOR  PTCI  FAILURE  TIMES 

',F10.2 

7 

/IX, 'LAMBDA  USED  FOR  PTC2  FAILURE  TIMES 

',F10.2 

8 

/IX, 'MEAN  SERVICE  TIME  FOR  OFT1  REPAIR 

',F10.2 

9 

/IX, 'MEAN  SERVICE  TIME  FOR  OFT2  REPAIR 

•,F10.2 

A 

/IX, 'MEAN  SERVICE  TIME  FOR  WSTI  REPAIR 

',F10.2 

B 

/IX, 'MEAN  SERVICE  TIME  FOR  UST2  REPAIR 

',F10.2 

C 

/IX, 'MEAN  SERVICE  TIME  FOR  PTTI  REPAIR 

',F10.2 

D 

/IX, 'MEAN  SERVICE  TIME  FOR  PTT2  REPAIR 

',F10.2 

E 

/IX, 'MEAN  SERVICE  TINE  FOR  PTCI  REPAIR 

',F10.2 

F 

/IX, 'MEAN  SERVICE  TIME  FOR  PTC2  REPAIR 

',F10.2) 

WRITE(6.12)  SGM0F1.SGH0F2,SGMWS1,SGMWS2.SGMPTI,SGMPT2 
1SGMPC2,NMTPER 

SGMPC1, 

FORMAT (IX  'STANDARD  DEVIATION  FOR  OFT1  REPAIR 

',F10.2 

1 

/IX, 'STANDARD  DEVIATION  FOR  OFT2  REPAIR 

',F10.2 

2 

/IX, 'STANDARD  DEVIATION  FOR  WSTI  REPAIR 

',F10.2 

3 

/IX, 'STANDARD  DEVIATION  FOR  WST2  REPAIR 

',F10.2 

4 

/IX, 'STANDARD  DEVIATION  FOR  PTTI  REPAIR 

',F10.2 

5 

/IX, 'STANDARD  DEVIATION  FOR  PTT2  REPAIR 

',F10.2 

6 

/IX, 'STANDARD  DEVIATION  FOR  PTCI  REPAIR 

',F10.2 

7 

/1X, 'STANDARD  DEVIATION  FOR  PTC2  REPAIR 

',F10.2 

8 

/IX, 'NUMBER  OF  MAINTENANCE  PERSONNEL 

',110.0) 

WRITE(6.13)HQ,OPROYI,OPRDY2,OPRDY3,OPROY4,OPROT5,OPROY6, 

10PRDY7,0PRDY8,F,CL0CX,TNREP 

FORMATdX, 'MAXIMUM  QUEUE  LENGTH 

',110.0 

2 

/1X,'OFTl  OPERATIONAL  READY  RATE 

',F11.0,'X‘ 

3 

/1X,'OFT2  OPERATIONAL  READY  RATE 

',F11.0,'X' 

4 

/IX, 'WSTI  OPERATIONAL  READY  RATE 

',F11.0,'X' 

5 

/1X,'WST2  OPERATIONAL  READY  RATE 

',F11.0,'X' 

6 

/IX, 'PTTI  OPERATIONAL  READY  RATE 

',FU.0,'X' 

7 

/1X,'PTT2  OPERATIONAL  READY  RATE 

',FU.0,'X' 

8 

/IX, 'PTCI  OPERATIONAL  READY  RATE 

',FU.0,'X' 

9 

/1X,'PTC2  OPERATIONAL  READY  RATE 

',FU.0,'X' 

A 

/1X, 'NUMBER  OF  MAINTENANCE  ACTIONS  THAT  TOOK 

1 

$ 

B 

/IX, 'MORE  THAN  THREE  (3)  HOURS  TO  CORRECT 

',F10.2 

C 

/IX, 'SIMULATION  RUN  TIME 

',F10.2 

D  /IX, 'NUMBER  OF  MAINTENANCE  ACTIONS  DURING  SIMULATION 

'110.0) 

RETURN 

END 
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EXPONENTIAL  RANDOM  VARIATE  GENERATOR 

DEFINITION  OF  VARIABLES 

ALOG  -  Natural  log  function 

DSEED  -  Seed  for  random  U(0,1)  ntinber  generator 

EXPON  -  Exponential  random  variate 

FMEAN  -  Mean  parameter  of  exponential  distribution 

GGUBFS  -  U(0,1)  random  nuiRier  generator 

R  ■  Random  nuiiber 


FUNCTION  EXPON  (FMEAN) 

DOUBLE  PRECISION  DSEED 
DATA  DSEED  /43/ 

GENERATE  A  U(0,1)  RANDOM  NUMBER 
RsGGUBFS(DSEED) 

»**1ltlHi»f*M*»»t»»»»*»ii*****»1l»»*»»**1l*****»»***»»*»***1l»»1Hl< 

GENERATE  AN  EXPONENTIAL  RANDOM  VARIABLE  UITH  MEAN  FMEAN 


EXPON=-FHEAN*ALOG(R) 

RETURN 

END 

r*********##«**«****** 


NORMAL  RANDOM  VARIATE  GENERATOR 

DEFINITION  OF  VARIABLES 

MEAN  •  Mean  of  normal  distribution 

SIGMA  -  Standard  deviation  of  normal  distribution 
NSEED  -  Seed  for  U(0,1)  random  number  generator 
NORML  -  Normal  random  variate 

K  -  Index  to  determine  which  method  of  computation  to  use 

PI  -  Pi 

RONE  -  A  U(0,1)  random  number 

RTWO  -  Another  U<0,1)  random  nuiber 

ZONE  *  A  N(0,1}  random  variable 

ZTUO  -  Another  N(0,1)  random  variable 


FUNCTION  NORML  (MEAN, SIGMA) 

REAL  ME AN, SIGMA, NORML 

DOUBLE  PRECISION  NSEED 

DATA  K  /O/,  PI  /3, 14159/,  NSEED  /58/ 

•*«*«««««•««*«««*«««*****•*****«•*«•**«**********«•***■ 

CHECK  TO  SEE  WHICH  N(0,1)  RANDOM  VARIABLE  TO  USE 

******#*«**«******«*«•****•**********•*•*****•********' 


GENERATE  TWO  U(0,1)  RANDOM  NUMBERS 

RONE ^GGUBFS( NSEED) 
RTWO^GGUBFS(NSEED) 


GENERATE  TWO  N(0,1)  RANDOM  VARIABLES 

ZONE=SORT(-2»ALOG(RONE))*COS(2*PI»RTWO) 

ZTWO=SORT(-2*ALOG(ROWE))*SIN(2*PI*RTWO) 


COMPUTE  NORMAL  RANDOM  VARIATE  WITH  PARAMETERS  (MEAN, SIGMA) 
FOR  MEAN  AND  STANDARD  DEVIATION 

*********************************  A****************************** 

NORML*ZONE*SIGMA+MEAN 

K=l 

RETURN 
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COMPUTE  NORMAL  RANDOM  VARIATE  WITH  PARAMETERS  (MEAN , S I GMA»*2 ) 
FOR  MEAN  AND  VARIANCE 

10  NORML=ZTWO*SIGMA+MEAN 
K=0 
RETURN 
END 
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-  GGUBFS 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


IMSt.  RCXJTINE  NAME 


COMPUTER  -  IBM/SINGLE 

LATEST  REVISION  -  JUNE  1,  1980 

PURPOSE  -  BASIC  UNIFORM  (O.I)RANDOMNUMBER 

GENERATOR- FUNCTION  FORM  OF  GGUBS 

USAGE  -  FUNCTION  GGUBFS  (OSEED) 

ARGUMENTS  GGUBFS  -  RESULTANT  DEVIATE. 

DSEED  -  INPUT/OUTPUT  DOUBLE  PRECISION  VARIABLE 
ASSIGNED  AN  INTEGER  VALUE  IN  THE 
EXCLUSIVE  RANGE  (l.DO.  2U74S3647.D0 
DSEED  IS  REPLACED  BY  A  NEW  VALUE  TO 
USED  IN  A  SUBSEQUENT  CALL. 

PRECISION/HARDWARE  -  SINGLE/ALL 

REOD.  IMSL  ROUTINES  -  NONE  REQUIRED 

NOTATION  -  INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 
INTRODUCTION  OR  THROUGH  IMSL  ROUTINE 
UHELP 

COPYRIGHT  -  1978  BY  IMSL,  INC.  ALL  RIGHTS  RESERVED 

WARRANTY  -  IMSL  WARRANTS  ONLY  THAT  IMSL  TESTING 

HAS  BEEN  APPLIED  TO  THIS  CODE.  NO 
OTHER  WARRANTY,  EXPRESSED  OR  IMPLIED 
IS  APPLICABLE. 


REAL  FUNCTION  GGUBFS  (DSEED) 

SPECIFICATIONS  FOR  ARGUMENTS 
DOUBLE  PRECISION  DSEED 

SPECIFICATIONS  FOR  LOCAL  VARIABLES 
DOUBLE  PRECISION  02P31M,D2P31 
D2P31M={2**31)  -  1 

D2P31  =(2**31)(OR  AN  ADJUSTED  VALUE) 
DATA  02P31M/2147483647.DO/ 

DATA  D2P31  /2147483648.D0/ 

FIRST  EXECUTABLE  STATEMENT 
OSEED  =  DMOO(  16807. DO  DSEED, 02P31M) 

GGUBFS  =  OSEED  /  02P31 

RETURN 

END 


GGUN0010 

GGUN0020 

GGUN0030 

GGUN0040 

GGUN0050 

GGUN0060 

GGUN0070 

GGUN0080 

GGUN0090 

GGUN0100 

GGUN0110 

GGUN0120 

GGUN0130 

GGUN0140 

GGUN0150 

GGUN0160 

GGUN0170 

GGUN0180 

GGUN0190 

GGUN0200 

GGUN0210 

GGUN0220 

GGUN0230 

GGUN0240 

GGUN0250 

GGUN0260 

GGUN0270 

GGUN0280 

GGUN0290 

GGUN0300 

GGUN0310 

GGUN0320 

GGUN0330 

GGUN0340 

GGUN0350 

GGUN0360 

GGUN0370 

GGUN0380 

GGUN0390 

GGUN0400 

GGUN0410 

GGUN0420 

GGUN0430 

GGUN0440 

GGUN0450 

GGUN0460 

GGUN0470 

GGUN0480 

GGUN0490 

CGU00500 

GGUN0510 

GGUN0520 
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APPENDIX  B 


SIMULATION  MODEL  OF  AIRCRAFT  SIMULATOR  MAINTENANCE 


LAMBDA  USED  FOR  0FT1  FAILURE  TIMES  25.80 

LAMBDA  USED  FOR  OFT2  FAILURE  TIMES  25.80 

LAMBDA  USED  FOR  WST1  FAILURE  TIMES  22.64 

LAMBDA  USED  FOR  WST2  FAILURE  TIMES  22.64 

LAMBDA  USED  FOR  PTT1  FAILURE  TIMES  28.05 

LAMBDA  USED  FOR  PTT2  FAILURE  TIMES  28.05 

LAMBDA  USED  FOR  PTC1  FAILURE  TIMES  21.99 

LAMBDA  USED  FOR  PTC2  FAILURE  TIMES  21.99 

MEAN  SERVICE  TIME  FOR  OFT1  REPAIR  1.14 

MEAN  SERVICE  TIME  FOR  OFT2  REPAIR  2.08 

MEAN  SERVICE  TIME  FOR  UST1  REPAIR  1.46 

MEAN  SERVICE  TIME  FOR  UST2  REPAIR  1.46 

MEAN  SERVICE  TIME  FOR  PTTi  REPAIR  2.23 

MEAN  SERVICE  TIME  FOR  PTT2  REPAIR  2.23 

MEAN  SERVICE  TIME  FOR  PTCl  REPAIR  1.05 

MEAN  SERVICE  TIME  FOR  PTC2  REPAIR  1.05 

STANDARD  DEVIATION  FOR  OFTI  REPAIR  0.22 

STANDARD  DEVIATION  FOR  OFT2  REPAIR  0.36 

STANDARD  DEVIATION  FOR  UST1  REPAIR  0.43 

STANDARD  DEVIATION  FOR  WST2  REPAIR  0.43 

STANDARD  DEVIATION  FOR  PTTI  REPAIR  0.74 

STANDARD  DEVIATION  FOR  PTT2  REPAIR  0.74 

STANDARD  DEVIATION  FOR  PTCl  REPAIR  0.07 

STANDARD  DEVIATION  FOR  PTC2  REPAIR  0.07 

NUMBER  OF  MAINTENANCE  PERSONNEL  2 

MAXIMUM  QUEUE  LENGTH  1 

OFTI  OPERATIONAL  READY  RATE  97. X 

OFT2  OPERATIONAL  READY  RATE  87. X 

WST1  OPERATIONAL  READY  RATE  97.X 

WST2  OPERATIONAL  READY  RATE  96. X 

PTTI  OPERATIONAL  READY  RATE  95. X 

PTT2  OPERATIONAL  READY  RATE  93. X 

PTCl  OPERATIONAL  READY  RATE  98. X 

PTC2  OPERATIONAL  READY  PATE  98. X 

NUMBER  OF  MAINTENANCE  ACTIONS  THAT  TOOK 

MORE  THAN  FOUR  (4)  HOURS  TO  CORRECT  0.00 

SIMULATION  RUN  TIME  6000.00 

NUMBER  OF  MAINTENANCE  ACTIONS  DURING  SIMULATION  550 
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SIMULATION  MODEL  OF  AIRCRAFT  SIMULATOR  MAINTENANCE 


LAMBDA  USED  FOR  0FT1  FAILURE  TIMES  25. 

LAMBDA  USED  FOR  OFT2  FAILURE  TIMES  25. 

LAMBDA  USED  FOR  WST1  FAILURE  TIMES  22. 

LAMBDA  USED  FOR  WST2  FAILURE  TIMES  22. 

LAMBDA  USED  FOR  PTT1  FAILURE  TIMES  28. 

LAMBDA  USED  FOR  PTT2  FAILURE  TIMES  28. 

LAMBDA  USED  FOR  PTC1  FAILURE  TIMES  21. 

LAMBDA  USED  FOR  PTC2  FAILURE  TIMES  21. 

MEAN  SERVICE  TIME  FOR  OFT1  REPAIR  1. 

MEAN  SERVICE  TIME  FOR  OFT2  REPAIR  2. 

MEAN  SERVICE  TIME  FOR  WST1  REPAIR  1. 

MEAN  SERVICE  TIME  FOR  WST2  REPAIR  1. 

MEAN  SERVICE  TIME  FOR  PTT1  REPAIR  2. 

MEAN  SERVICE  TIME  FOR  PTT2  REPAIR  2. 

MEAN  SERVICE  TIME  FOR  PTCI  REPAIR  1.05 

MEAN  SERVICE  TIME  FOR  PTC2  REPAIR  1.05 

STANDARD  DEVIATION  FOR  OFTI  REPAIR  0.22 

STANDARD  DEVIATION  FOR  OFT2  REPAIR  0.36 

STANDARD  DEVIATION  FOR  WST1  REPAIR  0.43 

STANDARD  DEVIATION  FOR  WST2  REPAIR  0.43 

STANDARD  DEVIATION  FOR  PTT1  REPAIR  0.74 

STANDARD  DEVIATION  FOR  PTT2  REPAIR  0.74 

STANDARD  DEVIATION  FOR  PTCI  REPAIR  0.07 

STANDARD  DEVIATION  FOR  PTC2  REPAIR  0.07 

NUMBER  OF  MAINTENANCE  PERSONNEL  3 

MAXIMUM  QUEUE  LENGTH  1 

OFTI  OPERATIONAL  READY  RATE  97. X 

OFr2  OPERATIONAL  READY  RATE  86. X 

WST1  OPERATIONAL  READY  RATE  96. X 

WST2  OPERATIONAL  READY  RATE  95. X 

PTT1  OPERATIONAL  READY  RATE  94. X 

PTT2  OPERATIONAL  READY  RATE  95. X 

PTCI  OPERATIONAL  READY  RATE  98. X 

PTC2  OPERATIONAL  READY  PATE  98. X 

NUMBER  OF  MAINTENANCE  ACTIONS  THAT  TOOK 

MORE  THAN  FOUR  (4)  HOURS  TO  CORRECT  0.00 

SIMULATION  RUN  TIME  6000.00 

NUMBER  OF  MAINTENANCE  ACTIONS  DURING  SIMULATION  553 
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SIMULATION  MODEL  OF  AIRCRAFT  SIMULATOR  MAINTENANCE 


LAMBDA  USED  FOR  0FT1  FAILURE  TIMES  25.80 

LAMBDA  USED  FOR  OFT2  FAILURE  TIMES  25.80 

LAMBDA  USED  FOR  WST1  FAILURE  TIMES  22.64 

LAMBDA  USED  FOR  WST2  FAILURE  TIMES  22.64 

LAMBDA  USED  FOR  PTT1  FAILURE  TIMES  28.05 

LAMBDA  USED  FOR  PTT2  FAILURE  TIMES  28.05 

LAMBDA  USED  FOR  PTC1  FAILURE  TIMES  21.99 

LAMBDA  USED  FOR  PTC2  FAILURE  TIMES  21.99 

MEAN  SERVICE  TIME  FOR  OFT1  REPAIR  1.14 

MEAN  SERVICE  TIME  FOR  OFT2  REPAIR  2.08 

MEAN  SERVICE  TIME  FOR  WST1  REPAIR  1.46 

MEAN  SERVICE  TIME  FOR  UST2  REPAIR  1.46 

MEAN  SERVICE  TIME  FOR  PTT1  REPAIR  2.23 

MEAN  SERVICE  TIME  FOR  PTT2  REPAIR  2.23 

MEAN  SERVICE  TIME  FOR  PTCI  REPAIR  1.05 

MEAN  SERVICE  TIME  FOR  PTC2  REPAIR  1.05 

STANDARD  DEVIATION  FOR  OFTI  REPAIR  0.22 

STANDARD  DEVIATION  FOR  OFT2  REPAIR  0.36 

STANDARD  DEVIATION  FOR  WST1  REPAIR  0.43 

STANDARD  DEVIATION  FOR  WST2  REPAIR  0.43 

STANDARD  DEVIATION  FOR  PTT1  REPAIR  0.74 

STANDARD  DEVIATION  FOR  PTT2  REPAIR  0.74 

STANDARD  DEVIATION  FOR  PTCI  REPAIR  0.07 

STANDARD  DEVIATION  FOR  PTC2  REPAIR  0.07 

NUMBER  OF  MAINTENANCE  PERSONNEL  4 

MAXIMUM  QUEUE  LENGTH  0 

OFT1  OPERATIONAL  READY  RATE  97. X 

OFT2  OPERATIONAL  READY  RATE  87.X 

WST1  OPERATIONAL  READY  RATE  98. X 

WST2  OPERATIONAL  READY  RATE  96. X 

PTT1  OPERATIONAL  READY  RATE  94. X 

PTT2  OPERATIONAL  READY  RATE  94. X 

PTCI  OPERATIONAL  READY  RATE  98. X 

PTC2  OPERATIONAL  READY  RATE  98. X 

NUMBER  OF  MAINTENANCE  ACTIONS  THAT  TOOK 

MORE  THAN  FOUR  (4)  HOURS  TO  CORRECT  0.00 

SIMULATION  RUN  TIME  6000.00 

NUMBER  OF  MAINTENANCE  ACTIONS  DURING  SIMULATION  552 
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